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CONTINUITY AND DISCONTINUITY OF THE BOUNDARY LAYER TAIL 


WILLIAM M. FELDMAN, INWON C. KIM 


Abstract. We investigate the continuity properties of the homogenized boundary data g for oscillating 
Dirichlet boundary data problems. The homogenized boundary condition arises as the boundary layer 
tail of a problem set in a half-space. The continuity properties of this boundary layer tail depending on 
the normal direction of the half space play an important role in the homogenization process in general 
bounded domains. We show that, for a generic non-rotation-invariant operator and boundary data, g 
is discontinuous at every rational direction. In particular this implies that the continuity condition of 
Choi and Kim [16] is essentially sharp. On the other hand, when the condition of [16] holds, we show 
a Holder modulus of continuity for g. When the operator is linear we show that g is Holder-^ up to a 
logarithmic factor. The proofs are based on a new geometric observation on the limiting behavior of g at 
rational directions, reducing to a class of two dimensional problems for projections of the homogenized 
operator. 


1. Introduction 


To motivate the questions considered in this paper, let us start by discussing the homogenization of 
oscillating Dirichlet boundary data problems, 

f F(D 2 u e ,^) = 0 in U 

\ ( 1 . 1 ) 

[ u £ = g( f) on dU. 

Here F(M,y ) is uniformly elliptic and positively 1-homogeneous in M, g{y) is continuous, and both are 
Z d -periodic in y. In the linear case we are considering operators of the form F(M, y) = —Tr (A(y)M) 
with 1 < A(y) < A and Z d -periodic in y. There is no problem to include a large scale x dependence in 
g but we omit it here for clarity. 


This type of problem has a singular behavior near boundary points x with inward normal direction 
v x aligned with a Z d -lattice vector, called rational directions. In order to mitigate the effects of the 
singularities the bounded domain U C is typically assumed to be uniformly convex, although more 
general assumptions which rule out large flat portions of dU are also be sufficient for the results discussed 
below, see [16,19]. In such domains it is known due to Feldman [19] that there exists g : —► R, 

continuous at irrational directions, so that u £ converges to u locally uniformly in U where u is the unique 
solution of, 

[ F(D 2 u) = 0 in U 

{ ( 1 . 2 ) 

\ u = g(v x ) on dU, 

where, again, v x is the inward normal of U at x G dU. Similar results have been obtained for linear 
divergence form equations starting with the work of Gerard-Varet and Masmoudi [21,22] and continued 
by several authors [2-4,28] . 

In this paper we study the continuity properties of the homogenized boundary data g by investigating 
the associated cell problem (1.4). Besides being a natural question on its own, continuity properties 
of g play an important role in obtaining rates of convergence for (1.2). In fact, if we could obtain 
Lipschitz continuity of g for the linear problem then we could also obtain an optimal rate of convergence 
that matches the rate for Laplacian operator. To our knowledge this particular connection has not 


2010 Mathematics Subject Classification. 35J60, 35J57, 35B27, 76F40. 

Key words and phrases, homogenization, oscillating boundary data, fully nonlinear elliptic equations, boundary layers. 


1 



WILLIAM M. FELDMAN, INWON C. KIM 


been written down explicitly in the literature, however it is implicit in the methods used in several 
works [20-22]. Let us point out that the typical strategy to study homogenization of (1.1) is by ensuring 
that the impact coming from singular boundary points are negligible: this is because in the linear case 
zero measure sets are not seen by the Poisson kernel, and in the nonlinear case an analogous argument 
applies for boundary sets of small Hausdorff dimension. In contrast, here we investigate the behavior 
of g as v x approaches rational directions. In the linear case we show, interestingly, that g extends 
continuously to the rational directions, and in the nonlinear case we show that discontinuity is generic. 

In the Neumann case the continuity of the corresponding g has been studied by Choi-Kim-Lee and 
Choi-Kim [15,16]. There it was shown that when the averaged operator F is rotation invariant, homog¬ 
enization holds and the homogenized boundary data is continuous. Following these works [19] showed 
homogenization for general F in the Dirichlet setting, due to the new observation that (1.2) has a unique 
solution if the discontinuity set of g(v x ) on dU has sufficiently small Hausdorff dimension. This brings 
up the natural question of whether the homogenized boundary condition could in fact be discontinuous 
when F is not rotation-invariant. Our main results are (?) an explicit estimate on the mode of conti¬ 
nuity for g when F is rotation invariant or linear, (ii) when F is not rotation invariant or linear, g is 
‘generically’ discontinuous at every boundary point with rational normal direction (see Theorem 1.3 and 
Corollary 1.4). These results seem to be new even in the linear case. 

We expect our main results in this paper to hold with parallel proofs in both the Dirichlet and Neumann 
case. On the other hand we hope to keep our illustration simple so that our main ideas are presented 
clearly. For this reason we will only discuss the Dirichlet problem, even though our arguments build on 
the framework introduced for the Neumann problem in [16]. We leave the task of proving parallel results 
for the Neumann problem, including the general homogenization results in [19], for the future work. 

We proceed to give a more precise, but still informal, derivation of (1.4) from (1.2). We begin by 
reminding the reader of the derivation of the cell problem determining g. We consider a rescaling of the 
solution u e of (1.1) near a boundary point xq £ dU with unit inner normal v Xo , 

v e (y) = u e {x 0 + £y). 

The limit of v e (Rv Xo ) as R —> oo and e —> 0, if it exists, will be the homogenized boundary data g{v x0 ) as 
long as eR —► 0. The behavior of v e outside of the oscillating boundary layer is the quantity of interest. 
To proceed with the analysis we inspect the equation solved by the v £ , 


F(D 2 v s ,y + e 1 x 0 ) = 0 
v E = g(y + £~ 1 x 0 ) 


in £ 1 {U — x o) 
on £~ 1 {dU — Xq). 


(1.3) 


P v ~{yv> 0} 
dP„. 


(1.4) 


Since F and g are assumed to be Z d periodic in y, £~ 1 Xq can be replaced by r e = £~ 1 xq mod Z d . Note 
that along various subsequences t e could converge to any r £ [0, l) d . This motivates the definition of 
the cell problem. Let v £ S’ 1 " 1 , r £ [0, l) d and if) be a continuous Z d -periodic function and define 
v u , T (-', (i/>, F )) : P v —» R. to solve, 

f F{D 2 v„, T ,y + r) = 0 in 
[ v ,,, T = i){y + t) on 

It is not too difficult to see, at least formally, that, 

\v E (y) - v Vxo>Ts (y;g{x 0 , -))| -* 0 as £ -> 0. 

From this identification we can replace understanding g(v) with the easier problem of understanding the 
limit v v<T {Rv) as R —> oo for every r £ [0, l) d . 

For irrational directions v the distribution of g on P v + tv is, in an appropriate sense, invariant with 
respect to r. For this reason it was possible to show, in [15,16,19], that, for irrational directions v, there 
exists a limit y(v, ip, F), the so-called boundary layer tail of v v ,ti such that 

sup sup \v v>T (y + Rv, ip) — ip, F)| —>• 0 as R —> oo. 

re[0,l)' i y<£dP v 


(1.5) 
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Note that in the context of (1.1) and (1.2) we should define g{v) := g{g,F,v). It was further shown 
in [15,16, 19] that g(y) : S d_1 \ RZ d — > R. is continuous, see Theorem 2.6. Thus the remaining question 
is to understand the limiting behavior of g{v) as v converges to a rational direction vq £ S ' d_1 D RZ d . 

The rate of convergence in (1.5) degenerates at rational directions where, in fact, the boundary layer 
tail does depend on r. Indeed the rational directions v are the possible discontinuity points of g. It 
turns out that the asymptotic behavior near the rational directions is actually quite structured and a 
more careful analysis is warranted. We will show that there is a multi-scale homogenization occurring, 
near-boundary in the micro-scale and then further away from the boundary in an intermediate scale, 
as irrational directions approach a rational direction. This phenomenon, partially described previously 
in [16], leads to a secondary homogenization problem associated with (1.4) with its own ‘cell problem’ 
and ‘effective operator’. This is far from obvious, and it will indeed be the main observation of the paper. 
Let us attempt to give a heuristic derivation of the secondary homogenization problem. The reader may 
wish to skip to the statement of the main results as the following description is unavoidably somewhat 
technical. 

We begin with a lattice point £ £ Z d \ {0} and its associated unit direction £. We may assume that 
£ is irreducible in the sense that the greatest common divisor of its entries gcd(£i, ...,£,*) is 1. A Z d - 
periodic ip on R d restricts to dP\ to be periodic with respect to a lattice on dPc with unit cell of size 
comparable to |£| (by the irreducibility). The limit of v^o{y + Rt;) as R — > oo exists by the periodicity of 
the boundary data. Again we refer to this limit as the boundary layer tail of v^ y Q. The same argument 
applies to v^ yT but unlike in the case of irrational boundary normal, generically, the boundary layer tail 
of V£ iT is not the same as that of v^ y o unless ( dP,\ + r) mod Z d = dP,\ mod Z d . We can concisely write 
down the set of limit points as, 

m^{t; {ip, F)) := ^lim v^ t c{R^', {ip, F)) which is a ||y-periodic function on R. (1.6) 

Now consider an irrational direction v which is very close to £. For a fixed yo G dP u the boundary dP v is 
close to dP ,t + {yo ■ £)£ on a very large region of size ~ \v — £| _1 centered at yo and so the respective cell 
problem solutions are also close in a smaller region (See Figure 1 in section 4.1). This observation leads 
to the conclusion that, far from the boundary, v v averages similarly to the following two-dimensional 
problem: 


| F{D 2 w iyri ) = 0 in ^ ^ 

\ W(, v = m ( (y ■ rj; {ip, F)) on <9P ? , 

where ry is the “approaching direction of v to £”, more precisely ?y is the unique unit vector rj _L £ so that 
the geodesic on the unit sphere, leaving £ at time t = 0 with velocity —rj, reaches v before time t = n. 

The limit of the homogenized profile /r(u, ip, F) as v -A £ can then be identified in terms of the approaching 
direction 77 . In other words we can show that for given ip, F, £ and v approaching £, there is a directional 
limit L = L^{rf), with g the approaching direction of v to £ as defined above, such that 

fj,{v, ip, F) — L^{g) —> 0 as v —\ |. (1.8) 

Precise statement of this result (stated quantitatively in Proposition 4.10) is the following: 

Theorem 1.1. Let /3 £ (0,1) and ip £ C 0, ^{T d ). For any £ £ Z d \ {0}, irreducible, there exists a 
function L^{-) = L^{-;ip,F) on unit vectors tangent to S d ~ 1 at ^ and a mode of continuity, W|^| i( g, such 
that the following holds: 

I n{v{t),ip,F) - L i {i/{ 0))| < W|£| i/3 (| v{t) - u(0)|) 
for any v : [0,1) —> S d ~ 1 a unit speed geodesic with v{ 0) = 
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The first half of the paper is spent to rigorously justify this derivation of the secondary cell problem and 
to obtain a quantitative estimate on the asymptotics of v u near rational directions. When the effective 
operators L£ are constant for every rational direction the quantitative estimates allow us to derive an 
explicit modulus of continuity for the homogenized boundary condition. From the previous arguments 
in [16] and [19] we know that Lg are constant, for instance, when F is either rotation invariant or linear. 

The characterization of the asymptotic behavior near rational directions described in (1.7) and (1.8) 
also opens the possibility of proving discontinuity of /i. One would just need to show that L j can be 
non-constant for some operator F, boundary condition if and lattice vector £ £ Z d \{0}. This simplicity 
turns out to be somewhat deceptive, as the situations where we can actually compute the boundary 
layer tail in (1.7) is when either the boundary data is trivial or the operator is linear, and L £ is constant 
in those cases. Another natural case to consider is when the operators are extremal, but then they are 
rotation invariant and thus L^ is constant, p is continuous. While it is difficult to come up with a specific 
example, it turns out that a better point is to show that a generic operator and boundary data will 
result in non-constant L j. In fact we are able to argue that if L ^ were to be constant, then we would be 
able to find many nearby F', if' with L £ non-constant. The perturbation of the operator is monotone 
and hence intrinsically nonlinear, and is designed to affect F in one direction 77 _L f while leaving another 
direction rf _L 77 , £ unaffected. The existence of 77 , rf mutually orthogonal and orthogonal to £ requires 
d > 3, and we are only able to achieve the desired perturbation of F when F = F is homogeneous. In 
fact, since the perturbations we make are quite explicit, besides showing that discontinuity is a generic 
phenomenon one can also generate specific examples of (if, F) where discontinuity of /i occurs. 

1.1. Main Results. The operators F(M,y ) discussed below will be positively 1-homogeneous, uni¬ 
formly elliptic with ellipticity ratio A, and Z d periodic in 77 . When F is linear we will write F(M,y ) = 
—Tv(A(y)M). If we say that F is spatially homogeneous we mean that F = F has no y dependence. 
For more details on these assumptions see Section 2.2. 

First we state our result about continuity. More details can be found in Section 5, for the improved 
estimate in the linear case see Section 7. 

Theorem 1.2. Let d > 2 and F such that (i) the homogenized operator F is rotation invariant or 
(■ ii ) F is linear. Then there exists a = a(A) £ (0,1) such that for any /3 £ (0,1), if £ C°’^(T d ) and 
v, v' £ S d ~ x \ RZ d we have 

I n(v,if,F) - p(v',if,F) | < C(d, A, / 3 )||^|| c .o,/ 3 (T d)|i 7 - F\ a0/d . 

In case (ii) we have additionally, 

\K V , 4 F ) - ^04 4 F)\ < C(d, A, || A|| c 5 (T d ) )|| 7 /;|| C 7 (T d ) |i7 - v'\ 1/d [l + (log ) 3 ]. 

For linear, divergence form systems a mode of continuity for p(-, if, F) is obtained by Gerard-Varet and 
Masmoudi [22] on the set of Diophantine irrational directions. Our result on the other hand is based 
on the mode of continuity near rational directions, and the modulus of continuity we obtain is uniform 
on the entire sphere. In the linear case it may be possible to combine these two results, but we do not 
pursue this here. 

Next we state our result about discontinuity. The statement is not completely precise, see Section 6 for 
the full details. 

Theorem 1.3. For d > 3, there is a residual set (in the Baire category sense) of continuous boundary 
conditions and spatially homogeneous nonlinear operators ( if,F ) such that p(-,if,F) does not extend 
continuously at any rational direction. 

The following question is left open. 

Open Problem. Does Theorem 1.3 hold when (i) F is taken to be inhomogeneous or (ii) d = 2? 



CONTINUITY AND DISCONTINUITY OF THE BOUNDARY LAYER TAIL 


5 


The argument used in the proof of Theorem 1.3 appears to be insufficient to address the above questions, 
however we do believe that the theorem holds in both cases (i) and (ii). 

The above results can be easily translated in terms of the original problem (1.2), since g(x ) = fj,(v x ,g, F ) 
for x £ dU with v x £ S d ~ x \ WL d . The details can be found in [19]. 

Corollary 1.4. Let U C be bounded uniformly convex domain, and let g be as given in (1.2). Then 
the following holds: 

(a) Suppose F is rotation invariant or linear, then the homogenized boundary data g(x) extends to 
be Holder continuous on dU, with mode of continuity in v x as given in Theorem 1.2. 

(b) Let d > 3. Then there is a residual set of g and F in (1.1), in the sense of Theorem 1.3 such 
that g, and hence u as well, are discontinuous at every x £ dU with v x a rational direction. 

1.2. Literature. There has been a surge of recent interest in the homogenization of oscillating boundary 
conditions, both in the linear divergence form and nonlinear non-divergence form settings. These works 
have much in common but there are some key differences which necessitate differing approaches. 

The problem is first addressed in the book of Benssoussan, Papanicolaou and Lions [12], which considers 
linear divergence form operators with co-normal oscillating Neumann boundary condition in general 
domains with no flat sides. The case of an oscillating Dirichlet boundary condition remained mostly 
open for quite a long time. For linear, divergence form systems recent progress began with the works 
of Gerard-Varet and Masmoudi [21, 22] where they show homogenization of the oscillating Dirichlet 
boundary condition problem with an explicit rate of convergence in L 2 (U). In that setting they show 
that the cell problem homogenizes at normal directions satisfying a Diophantine condition and that the 
rate of convergence to the boundary layer tail is better than polynomial. Continuing this investigation, in 
the direction of improved rates of convergence, are the works of Aleksanyan, Sjolin and Shahgholian [2-4], 
They identify the expected optimal L p convergence rate in general domains and obtain this rate under 
certain assumptions on the inhomogeneity of the operator. In a slightly different direction is the work 
of Prange [28] which extends the results of [21, 22] to include all irrational directions. He shows that 
the convergence to the boundary layer tail can occur at an arbitrarily slow polynomial rate without the 
Diophantine assumption. Perhaps the most relevant work to our paper is a recent result of Aleksanyan [1] 
on the continuity of the homogenized boundary condition. He shows for layered media, where the 
operator is independent of translations in the direction, that the homogenized boundary condition is 
as regular as the boundary data ip away from a possible singular set on Xd = 0. Compared to his result, 
we do not rely on any structure assumption on the operator, but on the other hand we obtain only 
Holder-i continuity in the linear case. It should be remarked that our result is in the non-divergence 
setting, nonetheless it may be possible for our approach to carry over to the setting of linear systems. 

Next we discuss the nonlinear, non-divergence form operators. For nonlinear operators there are several 
significant differences from the linear case. Firstly, due to the blow up procedure leading to the cell 
problem, the operators in the cell problem will always be positively homogeneous and therefore non¬ 
smooth at 0 (or linear). This makes the cell problem inherently impossible to linearize and so no 
regularity estimates better than C 1,/3 (or C 2 ’^ in the convex case) should be expected. On the other 
hand, higher regularity seems to be essential to obtaining arbitrary polynomial rate of convergence to 
the boundary layer tail at irrational directions as was done in the linear case by [22]. For these reasons 
obtaining arbitrary polynomial rates of convergence for the cell problem seems quite difficult if not 
impossible in the nonlinear case. The second problem, explicated for the first time in this paper, is 
that the homogenized boundary condition can be discontinuous. For linear operators a discontinuous 
boundary condition does not pose such a serious issue because, by the Green’s function representation 
the interior values of the homogenized solution can be estimated by measure theoretic norms of the 
homogenized boundary condition. In the nonlinear case no such “boundary ABP” estimate is known 
and so the uniqueness and stability of the homogenized problem is at issue (see [ 9] for a partial resolution 
to this problem). In regards to the literature, most earlier works address the Neumann problem: some 
special cases were discussed in Arisawa [5] in the half-space setting with periodic boundary data, and 
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also by Tanaka [29] using probabilistic methods. More general results were proved later by Barles, Da 
Lio, Lions and Souganidis [10]. Only just recently the full problem in general domains was considered 
by Choi, Kim [16] and Choi, Kim and Lee [15], wherein they show continuity of the homogenized 
Neumann boundary condition for rotationally invariant operators. For the Dirichlet problem Barles and 
Mironescu [11] obtained homogenization in half-spaces for a general class of nonlinear operators. The full 
problem in general domains was then considered by Feldman in [19]. The random case was considered 
in Feldman, Kim and Souganidis [20]. We also mention the recent work by Guillen and Schwab [23], 
where the half-space Neumann problem has been formulated as an interior homogenization problem for 
nonlinear non-local operators. 

1.3. Outline of the Paper. In Section 2 we start with notations and preliminary results to be used later 
in the paper. In Section 3 we prove the exponential rate of convergence for the half-space cell problem, 
when the boundary data is periodic on the boundary. While the proof is relatively straightforward, 
our result appears to be new for nonlinear operators. In Section 4 we investigate the behavior of the 
homogenized boundary condition ip,F) as v approaches a rational direction £ with £ £ Z d \{0}. We 
derive a second boundary homogenization problem that governs the directional limits as v approaches 
£. In Section 5 we prove Theorem 1.2 as a consequence of estimates in Section 4 and the Dirichlet’s 
Theorem (Theorem 2.11). In Section 6 we show that, when F is nonlinear, /r is generically discontinuous 
(Theorem 1.3). Finally in Section 7 we show that when F is linear and if} is sufficiently regular /r(-, F) 
is Holder-^ continuous up to logarithmic factors. In the Appendix we prove an extension of the result 
in section 3, which we make use of in section 7. 

1.4. Acknowledgments. We would like to thank Charlie Smart and Jason Murphy for helpful dis¬ 
cussions. We would also like to thank the anonymous referees for their very detailed and thoughtful 
comments which have helped very much to improve the presentation of the paper. Finally we would like 
to thank the hospitality of the Institut Mittag-Leffler where part of this research was conducted. Both 
authors are supproted in part by NSF grant DMS-1300445. 

2. Preliminaries 

This section contains notational conventions, fixing of the assumptions on the pde operators, statements 
of previously known results, and proofs of several technical Lemmas. The material here will be used 
throughout the paper and we suggest that the reader refer back as needed to this section rather than 
begin a careful reading here. 

2.1. Notation. We denote the half space with inner normal v by P v = {y : y ■ v > 0}. For a vector 
e £ R d \ {0}, e is the unit vector in the same direction e = |e| _1 e. We will occasionally need to project 
a vector e onto the orthogonal complement of another vector / £ R d , this we denote, 

Hf± e = e- (e • /)/■ 

We say that a constant C > 0 is universal if it depends only on the cllipticity ratio A and the dimension 
d. These constants may change from line to line without comment. If we need to refer to a specific 
universal constant which is not changing between lines we may call it Cg or Ci. For two quantities A, B 
we write A < B if A < CB for a universal constant C. If C additionally depends on a parameter b 
which is not universal then we will write A B. 

We will work with the function spaces of Holder continuous functions C k,/3 (X) for k £ N U {0} and 
/3 £ (0,1] with ( X , d ) a complete separable metric space. Most often A' = T" = R" mod Z” with metric 
inherited from Euclidean distance on M rl . We will repeatedly use the Holder semi-norm and norm for 
(3 £ (0,1], for a (f> : X —> R, 

\<t>\CW(X) : = sup Y|— and Mco-f>(X) = sup 

x^y(z:X CL\X , y) H x£:X 
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On R n (or T n ) the norms for the higher order Holder spaces are defined inductively for k > 1 and 
P S (0,1] by, 

= ll < / , llc fc - 1 .i(R") + \D k (f)\ c o.f>(^ n y 

2.2. Uniformly elliptic operators and viscosity solutions. We will work in the class of fully non¬ 
linear uniformly elliptic equations. Let A4dxd be the class of d x d real symmetric matrices. For 
F : Mdxd ->Rwe say F is uniformly elliptic if there exist 0 < A < A so that, 

XTr(N) < F(M) - F(M + N) < ATr (N) for all M, N £ M dxd with N > 0. (2.1) 

Lastly we define the class of uniformly elliptic operators, 

S\,a = {F : M dxd -t R : (2.1) holds } 

We will assume the following on F. 

(i) There is some A > 0 so that F(-,y) £ «Si,a for all y £ T d . 

(ii) F is Lipschitz continuous in y, 

\F(M, y ) - F(M, z)\ < <7(1 + \\M\\)\y - z\. 

(iii) F is positively 1-homogeneous, 

F(tM,y) = tF(M,y ) for all t > 0. 

We note that under the above assumptions F(M , y) is in fact an Isaacs operator arising from differential 
games (see for instance [13]), 

F(M,y) = inf sup — Tr(A ab (y)M) with 1 < A ab (y) < A. 

“e - 4 bes 

Next we recall the Pucci extremal operators associated with the ellipticity class S\, a, whose basic 
properties can be found in the book [13]: 

V+ a (M) := AE ei>0 ei + AX ei< oe* and P A A (M) := AE ei> o ei + AE ei< oej. 

Here e^’s denote the eigenvalues of M. The Pucci operators govern the worst possible behavior for 
viscosity solutions of F(D 2 i t, y) = 0 with F £ S\. a- More precisely note that for any M, N £ M dxd and 
any x £ R n we have 

— V + (M — N) < F(M, y) - F(N,y) < —V~(M — N). (2.2) 

It is not too difficult to check that the weak maximum principle holds for uniformly elliptic equations, 
the more difficult thing is the comparison principle. The following lemma, proved based the method 
of sup and inf-convolutions originally used by Jensen [24], shows that for uniformly elliptic nonlinear 
equations comparison principle for F follows from maximum principle for the Pucci operators. 


Lemma 2.1. Let LI be a domain in K n , and let F £ <S(A, A) satisfy (i). Let u and v satisfy, in the 
viscosity sense, 

F(D 2 u,y) < F(D 2 v,y) in LI. 

Then w = u — v satisfies, in the viscosity sense, —V AA (D 2 w) < 0 and —Vf A {D 2 w) >0 in LI. 

In addition to the role they play in the above lemma, the Pucci operators are useful because regularity 
results hold uniformly in the ellipticity class <S(A, A). The following result is from the book of Caffarelli 
and Cabre [13]: 


Lemma 2.2. Let F £ <S(A, A), and let u be a continuous viscosity solution of 

—V + {D 2 u ) < 0 and — V~(D 2 u) > 0 in B r { 0). 
Then for every a £ (0,1) there exists C = C(X, A, n, a) > 0 such that 


sup 

x,y£B T/2 (x 0 ) 


I u(x) -u(y) | < ^JL 


\x-y\ c 


sup u{x). 

r x£B r (x o) 
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Now given Lemma 2.1 we can discuss uniqueness/comparison principle in bounded domains and half 
spaces. To start let us consider a given bounded domain U C R ra with a smooth boundary and a 
continuous boundary data g : dU -A R. For an operator F which satisfies above assumptions (i)-(iii), 

f F(D 2 u, y) = 0 in U 

\ K ’ (2.3) 

\ u = g(y) on dU, 

We refer to [-3,17] for existence and uniqueness of viscosity solutions of (2.3), which is based on the 
following comparison principle. 

Lemma 2.3 (Comparison principle). Let u\, u 2 be viscosity sub- and supersolutions of (2.3) with 
boundary data g\ < g 2 . Then 

ui < u 2 in U. 

As for our cell problem (1.4) posed in the half-space, we will be using the following comparison principle 
for bounded viscosity solutions. 

Lemma 2.4 (Lemma 2.9, [19]). Suppose that U = P„ a half space with inward normal v £ S^ -1 . Let 
gi,g 2 £ C a (R n ) bounded and u\, u 2 be bounded sub and supersolutions of (2.3) with Dirichlet date g± 
and g 2 respectively, then 

ui < u 2 in U. 

A similar result will hold for sub/super solutions with sublinear growth, as one can easily check there is 
non-uniqueness once one allows for linear growth. 

2.3. Regularity in Two Dimensions. In d > 3 it is not known in general whether the solutions of 
fully nonlinear uniformly elliptic equations are smooth, examples of non-classical viscosity solutions in 
high dimensions (d > 12) have been given by Nadirashvili and Vladul; [25,26] . On the other hand in 
d = 2 it is a classical result of Nirenberg [27] that solutions are C 2,a for a small a. We will be able to use 
this result because the asymptotics near rational directions of the homogenized boundary condition in 
any dimension naturally turn out to be determined by a two-dimensional problem. We state the result 
using more modern terminology, but our statement follows easily from Nirenberg’s theorem in [27]. 

Theorem 2.5 (Nirenberg). There exists a(A) £ (0,1) and C( A) > 0 so that if u : Bi —> R is a viscosity 
solution of F(D 2 u) = f in Bi for some F £ and f £ C' 0 ” 8 (-B 1 ) then for a = min{ao, /3}, 

\\D 2 u\\co,«( B i/2 ) < C (A)[osc u + \\f\\co,S( Bl )}- 

£>1 

2.4. Results from Homogenization Theory. First we describe the results obtained in [19] regarding 
the cell problem, the Neumann counterpart is in [ 5,16]. Let lv, t (-; (if,F)) solve the cell problem (1.4). 
The following result says that, when v is irrational, v„ jT has a limit as y ■ v — ► 00 and the limit is 
independent of r. 

Theorem 2.6 (Theorem 1.2 of [19]). For v £ S' d_1 \ RZ d there exists F), called the boundary 

layer tail or homogenized boundary condition, such that, 

sup sup \v v , T (y + Rv) — g\ —> 0 as R —> 00 . 

re[o,i) d y£dP v 

Moreover p(-,if),F) is continuous on S d ~ x \RZ d . 

We will also need a rate of interior homogenization. In general this can be derived by the same methods 
used by Caffarelli-Souganidis [1 ] (also see Armstrong-Smart [6]). However in this paper we will only 
require an interior homogenization rate in the special situation where the solution of the homogenized 
problem in consideration is C 2 ’ a ° due to our two dimensional reduction and Theorem 2.5. In this case 
it is straightforward to obtain a rate of convergence, so we provide the proof. 
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For is £ S d 1 and R > 0, we consider the homogenization problem, 

( F(D 2 u e , f) = 0 in 0 < x v < R . ( F(D 2 u) =0 in 

< which homogenizes to < 

I u E = g{x) on x ■ v £ {0, R} I u = g(x) on 

where we are considering g : R d —>■ ffi. to be bounded and continuous. Suppose g satisfies 
g{x) = go(x ■ rj,x ■ v) for some unit vector rj L is and some go : M 2 —l ► K. 


0 < x ■ is < R 

x ■ is £ {0, i?}, 

(2.4) 

(2.5) 


Then by uniqueness u(x + t() = u(x) for any £ _L span {is,g}. In particular u(tis + sg) actually solves a 
fully nonlinear uniformly elliptic problem in d = 2 and hence has interior C 2 ’ a ° estimates by Theorem 2.5. 

We recall, for example from Evans [18], that for each M £ M.dxd there is a unique constant F(M) and 
a unique (modulo constants) Z d -periodic bounded solution v(y;M) of 

F(M + D 2 v,y) = F(M) in R d , (2.6) 

satisfying ||n(-; M)\\l°o < C(A,d)\\M\\. Again from [18], F turns out to be uniformly elliptic with the 
same ellipticity ratio A as F(M,y). 


Theorem 2.7. Let u E ,u : g be as given in (2.4) and (2.5). There exists 0 < a (A) < 1 such that for any 
/3 £ (0,1) and any R > 0, 

sup \u E {x) — u(x)\ < C(A, d)( osc g + R^\g\ c o,n){R~ 1 £) ald . 

0<x-v<R x-vG{0 ,R} 

Proof. After rescaling we may assume that R = 1 and U = {0 < x ■ v < 1}. Let S £ (0, ||s'||c°./ 3 (ai 7 )] t° 
be chosen later and u s solve 

( F(D 2 u s ) = <5 in U 

< J (2.7) 

[ u s = g{x) on dU. 

We claim that 

sup (u 15 — u\ < g<5. 

To prove this we look at w = u s —u which, by Lemma 2.1, is a solution of 

— V^~ a (D 2 w) < 5 < —Vf A (D 2 w) in U with w = 0 on dU. 

Comparing with 0 implies w > 0 and, for the other direction, let ip(x) = |(| — (x ■ is — |) 2 ), then 
— Vi A (D 2 ip) = 5 with ip > 0 on (?{7, so comparison principle implies w < ip < | in U. 

We will construct a supersolution barrier function based on u s to compare with u E away from the 
boundary. We begin by collecting uniform estimates on u s . Let 1 > h > 0 to be chosen small and call 
Uh = {x : d(x,dU) > h}. By the C 0 ’ 13 estimates up to the boundary - see Lemma 2.11 of Feldman or 
combine Lemma 2.2 above with Lemma 2.8 below - for both u and u E at unit scale, 

II«*IIcto.*( 1/) + ll ue ||c°>' 3 (U) < C'dls'll C°’t>(dU) + $) < C\\g\\co,P(du)i 
where we have used that S < ||9||c°>i 3 (a(7)- Moreover, due to Theorem 2.5 we have 
\D 2 u s (x)\ < C{h~ 2 osc u s + <5) < ChP- 2 \\g\\ c o,ft for x£Uh 

Bf da) 


where we have also used again S < ||<?||c°’£(d[/) anc l ^ 2 > 1. Similarly from Theorem 2.5, 

\D 2 u\x)\ c o, ao{Uh) < Ch^ 2+ ^\\g\\ c o.,. 


Note that for x £ U \ Uh there is y £ dU with \y — x\ < h and thus 

\u S (x) - u e (a:)| < |u 5 (a;) - g{y)\ + | u E (x) - g(y)\ < C\\g\\ C o.p {du) h p . 
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We wish to show that, in fact, the maximum of u e {x) — u s (x ) in U is obtained in U \Uh- Suppose 
otherwise, then there exists Xg £ Uh such that 

u e ( xo) — u S (x o) = ma x(u e — u s ). (2.8) 

In particular u 5 (x) + u e ( Xg) — u s (x o) touches u e from above at Xq. Let us define the barrier function 
4> e (x) := u e {x o) + Du s (x 0 ) ■ {x - x 0 ) + \(x - x 0 ) ■ D 2 u s (x 0 )(x - x 0 ) + e 2 u(§; D 2 u s (x 0 )) + ^\x - x 0 | 2 , 
where v is the corrector given in (2.6). Note that 

\4> e ( x o) - u e (x 0 )\ < C 0 \\D 2 u 5 \\ Lao{Uh) £ 2 . (2.9) 


One can verify that, using the uniform ellipticity and the definition of the corrector, 

F{D 2 (/) e (x ), f) > F{D 2 u 5 {x q)) - 6 > 0. 


Let us choose 

r = mm[e^ (\\D 2 u s \\ L o a(Uh) /\D 2 u s \co,o, 0 (u h ))^ ,h}. 
We claim that, for 5 sufficiently small and Co from (2.9), 

<f(x) > u%x) + 2C 0 ||DV|| i =o ( ^ ) £ 2 on dB r (x 0 ). 


( 2 . 10 ) 


The comparison principle then would yield that the same inequality holds in B r (x o), yielding a contra¬ 
diction to (2.9). We now verify that S can be chosen so that (2.10) holds. Using (2.8) we have 

<t> e (x) > u e (x) + -^x-xol 2 - C\\D 2 u 5 \\ Laa(Uh) £ 2 - C\D 2 u s \ c o,<*o(u h )\x - x 0 \ 2+ao , 


we have chosen r above so that the last two terms are of the same size on dB r (x o). Thus, evaluating 
this on dB r (xo) we have 


</> £ (z) > u £ (x) + Ae 2+ “°(||-D U \\ L -(u h )/\D u lc°.“o {u h )) 2+a ° - Ci\\D u \\ L °°(u h )£ ■ 


Now suppose r < h and let us choose 

$ < C\\g\\co.p max{e^khP-( 2+ao \hP- 4 e 2 }, ( 2 . 11 ) 

then due to the regularity estimates on u s given above we arrive at (2.10). If r = h then <5 = Mh^~ 4 e 2 
to get the same contradiction. By a parallel argument we can show that the same choice of 5 will result 
in the minimum of u e (x ) — u s (x) occurring in U \Uh- 

Now we put together the bounds obtained above. Since the maximum and minimum of W — u e are 
achieved in U\Uh for S as above, 

sup \u — u e \ < sup \u — u & \ + sup lu* 5 — u e \ < CS + C'||g||(7o./3/i /3 . 

U u u\u h 


Using S as chosen in (2.11), 

sup |u — u e \ < C\\g\\ C o,e(£^h 0 - {2+ao) + hP~ A £ 2 + h?). 
u 


2q Q 

By choosing h = e( 2 + a o) 2 we arrive at 

sup | u — u e \ < CWgWco.a^ l ' 2 +° : o ) ' 2 ■ 
u 


□ 
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2.5. Continuity up to the Boundary. We will use the following result repeatedly in what follows. 
It is a fundamental technical tool used in estimating the difference between cell problem solutions in 
nearby half-spaces. The result addresses the continuity-up-to-the boundary for solutions of the Dirichlet 
problem, but it can be also viewed as a localization result. 

Lemma 2.8. Suppose that uj : [0, oo) —> [0, oo) is a modulus of continuity and u < w(l) satisfies, 

i -Vt^D 2 u)< 0 in B 1 Ci K 

\ u{x) < w(|x|) on d(Bi fl K ) 

where K is any set satisfying 0 £ dK and B\ fl K C B+. Then there is a modulus w depending on A, d 
and u) such that 

u(x) < C(A,d)u>(\x\) in BiCiK. 

If u}(r) = r 13 for some j3 £ (0,1) then Q(r) = C( A, d,/3)r^. 

For us K will either be the upper half space P ed or an intersection of two half-spaces (see Lemma 4.4). 
Because it is not obvious how to calculate ui for general u> we will work with Holder continuous boundary 
conditions throughout the paper so that we get explicit estimates. The generalizations to arbitrary ui 
present only notational difficulties. 

Proof. By rescaling, without loss w(l) = 1. The proof is quite analogous to the standard barrier method 
for boundary continuity for harmonic functions, we just need to work with the Pucci operator instead 
of the Laplacian. Let </> be a positive, smooth function in B+ satisfying 

[ -Pf AD 2 4>) > 0 in Bi n K 

< ’ and (j)(x)<C 0 (A 1 d)\x\ in BiCiK. 

\ 1 on dB\ Cl K 

For example one can choose </> to be a rescaled translation of the downward pointing fundamental solution 
for the Pucci operator, 

4>(x) = L(1 - ]x + e d | 1-A(d-1) ) with L = ( min |1 - \x + e d | A(d " x) |)” 1 , 

]x. = l ,:c,j>0 

which one can check is actually a smooth solution of —P+ A (D 2 (j)) = 0 except at x = — e^. For each r > 0 
and an M > 1 to be chosen large consider the barrier, 

4> r {x) = uj(Mr) + ( sup u) + (j)((Mr)~ 1 x). 

Then <f> r (x ) > (sup BMrnK - u ) > u on dBMr H K and 

> <x(Mr) > w(|a;|) > u(x) on dl\ n Bm t , 

since ui is monotone. By comparison principle u < (f> r in Bmt Fl K and therefore it follows that 

( sup w(a ;))+<w(Mr) + ( sup u)+ sup < u(Mr) + ( sup u)+%, (2.12) 

B r DK B Mr nK x&B r r\K B Mr r\K 

where we have used that 4>(x) < Co|x| for the second inequality. To get a modulus of continuity, let 
e > 0, choose M > 2s~ 1 Cq and then choose r sufficiently small to make the right hand side in (2.12) 
less than or equal to e. 

On the other hand, since the argument is valid for every r > 0, applying the estimate repeatedly up 
until M n+1 r > 1, 

n— 1 

sup u{x) < 5>(A Pr)^) 2 - 1 + ( sup uUi^r- 1 
B r CK B\C\K 
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In case w(r) = for some /3 G (0,1), choose M 1 ^ = 2C'o(A, d) so that 

n—1 

sup u(®) < Mr 0 V + ( sup w) + (M-( 1 -^C' 0 )"“ 1 M-( n ~ 1 ^ 

B r nK B\f]K 

< 2Mr p + ( sup u)+2" (n - 1) M 2/3 r /3 

< C"(A,d,/?)(l + ( sup u)+)rA 

Biflif 


□ 


2.6. Some Number Theory. Lastly we present some elementary number theoretic results which we 
will make use of. When v G WZ d is a rational direction then there is some minimal T = T(y ) > 0 such 
that, 

(dP„ + Ti/) mod Z d = dP v . (2.13) 

Lemma 2.9. If £ £ Z d \ {0} is irreducible in the sense that gcd(£i,... , £<j) = 1 then, 

ni) = ir 1 - 

Proof. From Bezout’s identity there exists x G Z d so that £ • x = gcd(£) = 1. Then x ■ £ = |£| _1 and so 
x G dPs + ||y|. Thus, 

0 G [ OP £ mod Z d ] n [(9P| + i) mod Z% 

and so the two sets are the same. This shows that T(£) < A, for the other direction one just needs to 
note that £ • x is an integer for every x G Z d so that if £ • x 0 then |£ • x\ > 1. 

□ 

Lemma 2.10. If £ G Z d \ {0} then dP.'c is spanned by d — 1 vectors / J G Z d with |/ J "| < |£|. 

Proof. Without loss assume that |^| = arg m ax |£, ; | > 0 since £ ^ 0. Then call, for 1 < j < d — 1, 

l<i<d 


f 3 = — £,j£d and from the definition / J • £ = 0. 


□ 


Next we state a classical number theoretic result, the simultaneous version of Dirichlet’s approximation 
Theorem. The proof is by pigeon-hole principle. 

Theorem 2.11. For given real numbers an,..., a n and JVgH, there are integers pi, ...,p n , q G Z with 
1 < q < N such that 

1 

3. Asymptotics of Half-space Solutions with Periodic Boundary Conditions 

Here we consider the convergence rate for homogenization of half-space problems. First consider the 
solution v of the following problem in a half-space, 


(3.1) 


f F(D 2 v, y) = f{y) in P ed 
\ v = <f>(y) on dP ed . 

We assume that P, / and <j> are periodic with respect to linearly independent translations l\, ...Id-i G 
dP ed . We define the lattice of periodicity and its unit cell, 

2 := { z:z = J2 k :> G Z }’ : = | Vi : A j e [0,1) j, 


and L := diam(Q). 
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In this section we will only consider the case / = 0 for the simplicity of presentation. The proof of 
Lemma 3.1 for the general case / ^ 0, which is needed in Section 7, is presented in Appendix A. The 
calculations are rather delicate, since for later usage it will be important for us to keep track of the 
dependence on the unit cell size L. 

The following lemma states that the rate of convergence to the homogenized boundary condition will be 
exponentially fast depending on L and universal constants. This result is originally due to Tartar [30] for 
linear divergence form operators. To the best of our knowledge the result is new for nonlinear operators. 
The proof is an iterative argument using the ^-periodicity of the solution and the interior oscillation 
decay from Harnack inequality. 

Lemma 3.1. There exist n(<j>, F) and Co(A, d) > 0 such that, 

sup | v(y) — n\ < C( A, d)(osc(/)) exp(—coL~ 1 R), 

y-e d >R 

and this estimate gives the optimal rate up to the determination of cq ■ 


Proof. By rescaling we may assume without loss that osc <j> = 1. Let a(d, A) £ (0,1) be the Holder 
continuity exponent and Co(A,d) the constant in the interior Holder estimate for the maximal class 
(see [13]). For r > C]J a L we claim that 

osc v < C'o(L/r) afe for any k £ N. (3.2) 

Pe d +kred 

Supposing that this result holds, let k = [ R/r} to obtain 


osc v < osc v < exp I k log 


C 0 L c 


P c J +Red 


When we choose r = eC\j a L the estimate becomes, 

R 


osc v < exp — 


eC, 


l/c 


< C exp (— cR/L ) 


-l n~ x ! c 


with C = e and c = e i C' t 

It remains to prove (3.2) by induction. For k = 0 (3.2) follows from the maximum principle. Assuming 
(3.2) for k, we prove it for k + 1. Note that 

Vk(y) = Cg k r ak L~ ak v(y + kre d ) 

satisfies F(-, y + kre d ) = 0 in P ed with boundary data 4>k(y) = C^ k r ak L~ ak v(y+ke d ). Both the operator 
and the boundary data are periodic with respect to (£j)jZj translations by uniqueness, and osccfk < 1 
by the inductive hypothesis. Then by the interior Holder estimate in B r {re d ), 

Kl c»(B r , 2 (re d )) < C 0 r~ a and so osc v k < C 0 (L/r ) a , 

' Q+red 

where we have used r > 2L so that Q C H r / 2 (0). On the other hand v k is periodic on dP ed + re d with 
respect to the translations (Vj)jZi an d periodicity cell Q. Therefore 

osc v k < osc v k < C 0 (L/r ) a , 

Pc d +re d dPc d +re d 

where we have again used maximum principle for the first inequality. Rewriting this in terms of v, 

osc v = C£(L/r) ak osc v k < C k+1 (L/r) {k+1)a . 

Pe d +(k+l)re d Pc d +re d 

This completes the inductive proof. 

Lastly to show that the rate is optimal we take F = — A and = cos Then v(y) can be explicitly 
computed using separation of variables as 

v(y) = cos(^p-) exp(-^y 2 ). 
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Plugging in yi = 0 and 1/2 = R completes the proof since evidently p = 0 in this case. □ 

We will need a slight variant of Lemma 3.1 when the operator does not share the periodicity cell of the 
boundary data but its oscillations are at a smaller scale (see the proof of Proposition 4.3). We no longer 
assume that F shares the periodicity lattice of (f> and instead we suppose that there is 0 < e < L such 
that, 

for every y £ dP ed there is y' with | y — y'\ < e and F(M, ■ + y') = F(M , •) in P ed . (3.3) 

Lemma 3.2. Let <f>,v and L as given in Lemma 3.1 and F satisfies (3.3). Then there exists C,c> 0 
depending only on A, d such that, 

osc v(Re d + y) < C[(osc</>) exp(—cL _1 i?) + w„(e)]. 
y&dPe d 

Here lo v is the modulus of continuity of v at points of dP ed , 

u v (r):=sup{\v(y)-v(y')\:y£dP ed ,y'£P ed and \y-y'\<r}. 

Note that by Lemma 2.8 we have, for any fi £ (0,1), co v (e) < C(A, d, . 

Proof. The proof is a minor modification of the proof of Lemma 3.1 and so we mainly focus on the 
difference in the proof. We will prove that there is a(A, d) £ (0,1), Co (A, d) > 0 such that for r > 
(2C 0 ) 1 /«i 

osc v < Cq (L/r) ak (osc (f) + C 0 uj v (e) for any k £ N. (3.4) 

P Sd +kre d 

Following the proof of Lemma 3.1 we are done as long as we can show (3.4). 

The proof of (3.4) is again by induction. We wish to show that (3.4) holds for all k. Assuming that 
(3.4) holds up to k we prove for k + 1. Note that Vk(y ) = v(y + kre d ) solves the equation 

F k (D 2 v k , y) := F(D 2 v k ,y + kre d ) = 0 in P ed 

with boundary data 4> k (y) = v k -i(y + re d )• Note that F k satisfies the same assumption as F. By the 
interior oscillation decay for the ellipticity class, there exists a universal constant C\ > 1 such that 

osc v k < CAL/r) a osc (h k for any £ £ Z. 

C+Q+re d dP Cd 

For an arbitrary y £ dP ed let £ £ Z such that y £ £ + Q, 

\v k (y + re d ) - v k (re d )\ < 2 osc v k + \v k {£ + re d ) - v k (re d )\ 

z+Q+red 

< 2C' 0 (L/r) a osc fi k + \v k (£ + re d ) - v k {re d )\. (3.5) 

ap *d 

The second term on the right hand side above appears since v k is no longer periodic. By the assumption 
on F there is £' £ P ed with \£f — £\ < e and F(M , ■+£') = F(M , •) in P ed for all symmetric matrices M. 
Therefore we can estimate, 

| v k (£ + re d ) -v k (re d )\ < \v k (£ + re d ) - v k {£' + re d )\ + \v k (£' + re d ) - v k (re d )\. 

The first term can be estimated by, 

| v k {£ + re d ) - v k {£' + re d )\ < C 0 {\£ - £'\/r) a osc f> k . (3.6) 

dP -d 

For the second term it suffices to bound w(y) = v(y + £') — v(y) in P ed . From the invariance of the 
operator F under translation by £' and Lemma 2.1, w is a viscosity solution of 

-P+ a (D 2 w) < 0 < -VfJD 2 w) in P ed . 

Using the periodicity of cf> we have w(y) = v(y + £') — <f>(y + £) on dP ed . Then by maximum principle 
and the definition of uj v , 

sup |u>| < sup \w\ = sup | v(y + £') — c/)(y + £)\ < uo v {\£ — £'\). 

Pe, dPe, dPe, 


(3.7) 
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Plugging (3.6) and (3.7) into (3.5) and using \l — l' \ <£< L we obtain, 

\vk(y + re d ) -v k {re d )\ < 2,C' 0 {L/r) a osc 4> k + u v (e). 

dP - d 

We can conclude now since, 

osc v = osc v k 

dPe d + (k+l)re d dP ed +re d 

< 3Ci{L/r) a osc 4> k + w^(e) 

dP *i 

< 2,C l C^{L/ri k+l)a osc <f>+ (3Ci C$ {L/r) ka + l)w„(e) 

< C k+1 (L/rY k+1 ^ a osctp + Coco v (e) 

where the last step holds if we choose Co = 3Ci and r > 2 1 t° l LC^ a since 1 < C\. 


(3.8) 


□ 


4. Asymptotics of half-space solutions near rational directions 

In this section we study asymptotic behavior of half-space solutions as the normal direction v approaches 
a rational direction, £ for some £ £ Z d \ {0}. Let us recall the solution v v ^ T of the cell problem (1.4) 
defined for a direction v £ S d ~ x , a t £ R d and a continuous Z Aperiodic ip by, 

f F(D 2 v„ t , y + t) = 0 in P v 

{ (1.4) 

[ Vv, T = ip{y + t) on dP v . 

Due to Theorem 2.6, for irrational directions there exists a limit p(v,ip,F) such that, 

sup sup \v V}T (y + Ru; (ip, F)) — /x(u, ip, F)| — > 0 as R —> oo. (4.1) 

r6K d y£dP u 

We are interested to understand the asymptotic behavior of /i(u, F) as v approaches a rational direction 
£ for £ £ Z d \{0}. The limiting behavior, it will turn out, depends on the direction of tangential approach. 
The main result of this section (stated quantitatively in Proposition 4.10) is the following: 

Theorem 4.1. Let (3 £ (0,1) and 'll) £ C°’^(T d ). For any £ £ Z d \ {0}, irreducible, there exists a 
function L^(-) = L^(--,ip,F) on unit vectors tangent to S d ~ 1 at £ and a mode of continuity, u)£,p, such 
that the following holds: 

I h( p {t),if,F) - 0))| < W|£| )j8 (]i y(t) - u(0)|) 

for any v : [0,1) S d ~ x a unit speed geodesic with u(0) = 

The basic idea behind these asymptotics already appeared in [15, 16] for the Neumann problem as a 
part of their proof that, when F is rotation invariant, p(-,il>,F) has a continuous extension from the 
irrational directions to the entire unit sphere. The proof proceeds by a series of reductions which will 
be carried out by a multi-scale homogenization argument. Our analysis is a quantitative and improved 
version of the proof given in [16] in the following sense. First we have tried to obtain optimal estimates 
at each stage of the argument. We do this with the hope of clarifying the proof and of achieving 
improved quantitative results on the continuity of p, in the end. Secondly we introduce the directional 
limit L and observe that L depends on a two-dimensional projected version of the problem (see Section 
4.3). It is for this reason that we are able to use Nirenberg’s two dimensional regularity result and the 
corresponding interior homogenization result, Theorem 2.7. By this careful exposition we are able to 
obtain a precise characterization of the asymptotic behavior of /j at rational directions and its dependence 
on the operator F and boundary data ip. With this characterization we are able to understand both 
continuity and discontinuity, Sections 5 and 6 respectively, in a unified way. 
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4.1. Step 1: Replacing the Boundary Condition at an Intermediate Scale. Let CeZ d \ {0} 
be irreducible and let solve (1.4) in P ? . By the results of Section 2.6 the boundary data ip lap,, is 
periodic with respect to a lattice on dP^ with unit cell size < Cd\f\. The result of the previous section 
implies that for each r £ there is a limit at infinity in the £ direction. The limit for r, t' £ is the 
same when r, t' are both in dP^ + t£ modulo Z d for some i £ R. By Lemma 2.9 this is exactly when 
r • | mod ||yZ = t' ■ £ mod ||yZ. We define 

m 5 (t; (ip,F)) := Um^v^Ri) 

which is continuous, y|y-periodic on K. By definition we have 

lim v itT (Ri) = m ? (r • (■ ip,F )). 

it—^ oo 

Generally speaking inherits the up to the boundary regularity of the cell problem solutions. In the 
general fully nonlinear case this is limited to C 0,1 , but for linear operators the result would hold for 
arbitrary C' fc,/3 , see Section 7. 


Lemma 4.2. If ip is continuous with modulus u> then mj(-; ( ip,F )) is continuous with the new modulus 
Co from Lemma 2.8. In particular for j3 £ (0,1), 

ll TO ?(-; (i’,F))\\ C o.fi < C(A, d, f3)\\ip\\ C o,t>. 


Moreover we also have, 


||m 5 (-; {ip,F)) Hco.i < C(A,d,/3)\\ip\\ C i,0. 


The proof is a straightforward application of the boundary continuity estimates Lemma 2.8 combined 
with the definition of m^, and is postponed till the end of this section. We drop the dependence of mj 
on ( ip,F ) as long as there is no ambiguity. Due to Lemma 3.1 of the previous section, 

sup t |(-) — mc(f)| < C(oscip ) exp(—cP/|£|). 

P ( + Rp 

Let v £ S d_1 be an irrational direction and p,{v, ip, F) the boundary layer tail of the cell problem solutions 
Vv,t (see Theorem 2.6 for the definition of p). Since the limit is independent of r (from Theorem 2.6) we 
simply refer to v v = iv,o when v is irrational. We consider the asymptotics of ip, F ) as v approaches 

i 

When v ^ f there is a unique vector p !_ £ (see Figure 1) so that, 

v = (cos \r/\)f - (sin \ri\)f] with \r/\ = \u - || + 0(\v - £| 2 ) (4.2) 


It may be helpful, although it is not essential, to note that this is just the minus of the inverse of the 
exponential map exp^ : T^S d ~ 1 —► S ,d_1 \ {—|} where T/S^ -1 is the tangent space to S ,d_1 at £. The 
goal of this section is to show that, after moving to the interior and rescaling, the cell problem solution 
v v is very close, in terms of \v — £|, to solving, 


F (D 2 w i ^, \r)\ x y) = 0 in 

w t,v = m dv • v) ° n dP € , 


(4.3) 


in their common domain of definition. We do not claim that (4.3) has a boundary layer tail. Indeed 
the periodicity lattices of the boundary data and the operator may not be aligned. On the other hand, 
by Lemma 3.2, it will almost have a limit up to an error small in ( 77 ! and this will be sufficient for our 
purposes. More precisely we aim to prove: 


Proposition 4.3. Let f £ Z d \ {0}, irreducible, and v £ S d 1 \ RZ d with |£ — |£|i/| < 1/2 then, for any 

( 0 , 1 ), 

li{v,F,ip) - lhnmf w itV (Ri', Wh-F)) < C{A,d,P)\ip\ C o,f>\£ - log ■ 

The parallel statement holds for the lim sup as well. 



CONTINUITY AND DISCONTINUITY OF THE BOUNDARY LAYER TAIL 


17 


-11 




Figure 1. The appearance of m^ at an intermediate distance from dP v . 


We remark that the result does not depend on the Holder continuity of ip (any continuity modulus for 
ip would yield an analogous result). Furthermore when ip £ C 1 ’^ the estimate can be improved to 

ju(iy, F,ip) - liminf w itr/ (Ri; (ip,F)) < C(A, d, 0)\\ip\\ c i, fi |£ - |g|^| log |g _jg H ■ 

Let us give a heuristic proof of Proposition 4.3, which is illustrated in Figure 1. Pick a point yo £ dP v . 
In a neighborhood of yo the boundary data ip for v v is very close to that of ip restricted to dP,\ + (yo ■ £)£. 
This causes v v to be close to m^(yo ■ £) at Vi : = Vo + Rf f° r R = o(|£ — v \~ 1 )- Next, observe that 
yo ■ f ~ yo • rj, since f — y is almost v and y 0 is perpendicular to v. But now, since rj is perpendicular to £, 
we have yo ■ rj = yi ■ rj. Consequently one can now say v v (y) is now close to mc(y ■ rf) R -away from dP u , 
and this describes the near-boundary homogenization for v u . Now taking m^(y ■ rf) as the new boundary 
data for the interior homogenization, we arrive at the interior problem (4.3) and Proposition 4.3. 

The actual proof is slightly more involved for technical reasons. 

Lemma 4.4. Let f £ Z d \ {0} be irreducible. For |£ — |£|u| < 1/2 and /3 £ (0,1) let us define Rq := 

c“ 1 1 £ I log , ,, 1 -, . Then we have 
Iflk-fl 

sup \v v {Rov + y) — m^(y ■ £)| < C(A,d, P)\ip\ c o. P \f - |gM g log 

yedP u l?l 1 

We remark that the log term in above esimate can be improved slightly as may be noticed from the 
proof. 


Proof of Lemma 4-4- We first show that for all R > 1, 

sup | v v (Rv + y) -mt-fy-i) \ < (osc ip) exp(-oR/|£|) + \ip\ C o,pR p \v - ||' 3 . 

yGdPv 

This will imply the desired result by choosing R 0 = c~ 1 |^|log- ^ ^ and using (osc ip) < \ip\ C o,p. Fix 

yo £ dPv and we consider comparing v v with the solution w of, 

f F(D 2 w,y) = 0 in P^ + y 0 
1 w = ip{y) in dP^+yo- 

Note that yo £ dP / + (yo • f)f and therefore, using that v ■ f > 1/2, 

\w(y 0 + Rv) - m(y 0 ■ |)| < C(osc ip) exp (—c|^| _ 1 -R). 


(4.5) 
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On the other hand for y G d[(P^ + yo) fl P„] there exists y' G dP,\ + yo such that 

W -y\ < W - i\\y - yo\, 

and so by the Holder continuity of w up to the boundary, 

\w(y) - i>{y)\ < \w{y) - w{y')\ + \ ip(y') - iffy) I < C(d, N)\tp\ c o,p\v - £\P\y - yof. 

The same argument holds for v v and by combining the two estimates we have 

I v„(y) - w(y)| < min{C\’tp\ c o,e\v - £f\y - y 0 f ,osci/>} for y G <9[(P ? + y 0 ) O P v \, 
where the second term is from maximum principle, min^ < < ma xip. Now we claim that, 

My) ~w{y)\ < C(A,d,P)m.m{\ip\ C o,f>\v - £f\y - yo^ioscil)} for y G (P^ + yo) fl P v , (4.6) 
but this is just a rescaling of Lemma 2.8. In particular (4.6) combined with (4.5) implies, 

\v„(Rv + y 0 ) - m(y 0 ■ £)| < \v v (Ru + y 0 ) - w(Rv + y 0 )\ + C^lco^R 13 ^ - 

< (osc ip) exp(—cii/|£|) + \i/)\ 0 o,i>R fl \v - ||^. 

This was the desired estimate. □ 

Next we return to the proof of Proposition 4.3 from Lemma 4.4. 

Proof of Proposition 4.3. By maximum principle in the domain P„ + Rqv Lemma 4.4 implies that, 

My + Rqv) - u(y )I < - l|' d log j^jj^i for y G P v (4.7) 

where u solves, 

j F(D 2 u,y) = 0 in P„ 

[ u(y) = mt(y ■ i) = mz(y ■ U.„±£) on dP„, 

where we recall that := £— (£• v)v is the orthogonal projection onto dP u . In particular an estimate 

of the same form as (4.7) holds for the respective boundary layer tails. Call rj 0 = and recall that 
we had defined r) as 

77 := — expG 1 (^) defined so that v = (cos \y\)i — (sin \r)\)rj. 

From the definition of 77 we calculate, 

Vo = n v ±i = £-(£■ v)v= (sin \y\) 2 i + (sin|r?|)(cos \rj\)fj, 

and so, 

\Vo ~ v\ < I sin \v\ ~ \v\\ + \Vo - (sin|? 7 |)?)| 

< | sin |? 7 | - |r?|| + (sin \r]\)y/2(l - cos |r 7 |) 

<N 2 


Now we rescale to u(z) = u(|77| 1 z) which solves 

j F(D 2 u,\ri\- 1 z) = 0 in P u 
\ u(z) = m^(z ■ |?7|” 1 7?o) on 8P V , 

and estimate the difference of u and w^_ r/ in their common domain P v n P^. The aim is to obtain an 
estimate on the difference of their respective boundary layer tails. From here the proof will follow a 
familiar argument. From the estimates above and Lemma 4.2, 

I rnd z ■ \v\~ l Vo) - ™d z ’ v)\ < ■ 
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Using this we bound the difference u — u^. r/ for z £ d(P v D P$). First note that for 2 £ <9P„ fl P^ there 
is z' £ dP ,f with \z' — z\ = \z\\q\. Therefore we have 

Hz) - w i,ri(z)\ < | m^z • \vHvo) - mc(z ■ f})\ + \w i>ri (z) - m^z ■ fj)\ 

< \m^\ c o,p\z\ 0 \r]\^ + | w^z) - m^z' ■ fj)\ + \m^z' ■ fj) - m^z ■ fj)\ 

< Clm^co.olzfli^ 3 


where the middle term in the second line is estimated using the continuity up to the boundary of 
from Lemma 2.8. Combining this with osc< |mj| c o,/3and \mc\co,p < C(A, d)\ip\QO,p we have 

j —V^ k {D 2 {u — w^ v )) < 0 in P„nP 4 

\ {u-w (<v )(z) < C{A,d)\ip\ c o,fimm{\z\^\v- on d{P u nPe). 

Therefore by the rescaled version of Lemma 2.8, 

u(R£) — W£ tri (R£,) < C(d, A, f3)\ij)\co,pBP\v — ^|' 3 for any R > 0. (4.8) 

At this stage we want to combine this estimate with the exponential convergence of u , to their 
respective boundary layer tails, but there is a minor technical issue that F(M, \r/\~ 1 z) does not share 
the same periodicity lattice as m^{z ■ rf). However, the conditions of Lemma 3.2 do hold and the rate of 
convergence established in Lemma 3.2 combined with (4.8) implies, for any R > 0, 

limsupu(-R'£) — liminf w^ tV (R'^) < (oscto^) exp(—c|£|P) + \ip\ c o,n 

R'^oo R'-i-oo 

Now we are free to minimize over R > 0, then plugging in \r]\ < \v — £| to obtain 
limsup{t(P£) — liminf wc „(i?£) < C\ip\ c o,f> W ~ log 1 ?, ■ 

R-> 00 R->oo ’ l^-tl 


Finally combining with (4.7) and the remark below it that the same estimate holds for the boundary 
layer tails, 


v(v,il’,F) - liminf w^ v (R4) < C( A, d,ft) \ip\ c °>f>\tf W ~ i\ P log m l_£\ • 
A symmetric argument yields the same estimate for limsup^^^ u>£ iT) (P£) — ip, F). 


□ 


Proof of Lemma ^.D. For the purposes of this proof it will be useful to work with a slightly different 
definition of the cell problem solution. We call v^^ T (y) = uj jT (j/ — r) which now solves, 


F(D 2 v iiT ,y) = 0 in P ? + r 
V£ iT = ip(y) on dPc + t. 


(4.9) 


Of course the boundary layer tail remains unchanged. The point is that the now solve the same 
interior equation for all r £ but in different domains. When r — t' is small the domains are close and 
we can combine the boundary continuity estimate of Lemma 2.8 with comparison principle Lemma 2.4 
to estimate the difference of the cell problem solutions, and hence of their boundary layer tails as well. 
It suffices to estimate the continuity of at t = 0. Let /3 £ (0,1), e > 0 and any y £ dP^, 


~ k,o(v) = - Hy) 

= v£,- s i(v) - v t,-ei(y - £ i) + Hv - £ 0 - Hv) 

< e^(sup ||u{,r||c t) ./ 3 (P £ ) + IIV’llc°.' 5 (T‘i)) 

T 

< C(d, A, f))£^\\lp\\ C 0,ff(jdy 
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Then, by maximum principle the same estimate holds in P^ and therefore, 

I m^-e) - m ? (0)| = lim _ e f(P£) - 6f,o(P|)| 

it—>• OO s ’ s 

< sup |tL _AR£) - V(,o{R£)\ 

Pe 

<C(d,A,^M C o,, m - 

Parallel arguments work for e < 0. To get the Lipschitz estimate use the fact that 
sup||^, T ||coa ( p 5 ) < C(d, A, /3)\\'tp\\ C i,f>(jd' ) for any /? > 0. 

r 

□ 


4.2. Step 2: Interior Homogenization at the Intermediate Scale. From the reduction performed 
in the first step we are left to consider the following problem. For an rj _L £ with \y\ >0 small, 

f |) = 0 in P ? 

\ w^ v = m^(y ■ fj) on dP s , [ = m^y ■ y) on 

We wish to make this convergence quantitative so that we can get an estimate of the difference between 

/r(i/, F, ip) and boundary layer tail of the homogenized problem in (4.10). 


which homogenizes to 


F(D 2 w m) = 0 
w it f, = mt(y ■ fj) 


c)l\. 


(4.10) 


At this stage it is useful to note that the homogenized solution uJf.r) is actually two dimensional. The 
key observation here is that since the boundary data only varies in the fj direction and the homogenized 
operator is translation invariant, the solution w£ t fj only varies in the fj , £ directions. This is a simple 
consequence of uniqueness. 


Claim, only depends on x ■ £ and x ■ fj. 

We prove the claim only to emphasize the importance of passing from to w^y r 

Proof. For any ( _L fj, £ and t £ R note that w' = u>£,fj(y + K) solves 

F{D 2 W r ) =0 in P^ + £ = with w'(y) = m((y + tCf) • fj) = m{y • fj) on 9P ? . 
This is of course the same equation satisfied by so by the uniqueness of bounded solutions 


□ 


In particular we have reduced to a situation where, by Nirenberg’s Theorem, the homogenized solution is 
C 2 ’ a ° on the interior. By using the exponential rate of convergence to the boundary layer tail established 
in Section 3 combined with Theorem 2.7 we are able to show, up to a logarithmic factor, that the same 
rate of convergence holds for (4.10). 

Lemma 4.5. Let y 1 £ with |? 7 ||£| < 1/2. Then there is a(A) £ (0,1) such that for any /3 £ (0,1), 

\wt,v(y)- w t,v(y )I ^ C( A >d)|V>|c°’i 3 (Td)l£l^ (a ”* ) M a ^(log 

and, in particular, the same estimate holds between the boundary layer tail ofwand lim inf w^ v (Rf) 
or limsupp^ w^ v (R£). 

Before we proceed with the proof we state a consequence of Proposition 4.3 and Lemma 4.5. 

Lemma 4.6. Let £ £ T, d \ {0} irreducible and v an irrational direction with y = y{v) as in (4.2). Then 
there is a(A) £ (0,1) such that for any f3 £ (0,1), 

\n(v,ip,F) - lim w^ v {Ri)\ < C(A, d, P)\il>\c°.t>(T*) 
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Proof of Lemma 4-5- We would like to apply Theorem 2.7 to prove the Lemma, however a modified 
argument is necessary since and ^ are solutions in an entire half-space. In order to replace with a 
homogenization problem in a bounded width strip we use that converges with exponential rate 

to its boundary layer tail JL and u^ ;7) (/?,£;), although it does not quite have a boundary layer tail, converges 
with exponential rate to a neighborhood of width small in // centered at any of its subsequential limits 
/i (see Lemma 3.2). 

Now we begin with the technical details of the proof. First recall that is j|y-periodic on M. Therefore 
m z{y ' v) is y!y-periodic on dP^ in the direction fj and constant in the directions orthogonal to i). Due to 
Lemma 4.2 we can estimate 

osc mg < | £|P^ | TO£ | c'o.iS (®) <c\z\-e\i,\co.e, 

and 

l m £ lco>“/3(R) < C°’P( R)- 

Next let /.i and JL respectively denote any subsequential limit of and the limit of as 

R —► oo, and let a(A) be as given in Theorem 2.7. Then from Lemma 3.2 we have 

\ w £,ri(y) - v\ + \ w t,v(y) - vl <C(oscm^exp(-c\^\R) + C\m^\ c o, a p\r]\ a0 for y £ P i + R£. (4.11) 

We use (4.11) to restrict to a domain where we can use Theorem 2.7, then we simultaneously are able 
to estimate y — y and — w | i?) . Fix an R to be chosen and consider, 


u>£,r,(y) : = w a,v(y) + R VI 


(n-n)+ sup [k{,r,(-) - n\ + - Ml] 

dPz+Ri 


(4.12) 


Note that with this modification W£, v still solves the same equation as in P £ with the same boundary 
condition on dP^ but also 

^,r,(y)>w^ v (y) on dP^ + Rf. (4.13) 

Now Theorem 2.7 implies that 

wt,r,(y) - ^ + V)(r _ 1 m)V Vc-°./tr)- (4.14) 

Note that we are not quite applying Theorem 2.7 directly. To be precise we first solve the equation 
F(D 2 u , |^y) = 0 with boundary data matching w^ v in 0 < y ■ £ < R. By comparison principle and the 
ordering (4.13) we know u < u>£,n- On the other hand from Theorem 2.7 we have the desired estimate 
for |u — uJ£ ]77 |, combining these two steps we get (4.14). 

Rewriting (4.14) in terms of w^ v using (4.11), 

Wi,r, (y) - w^riiy) < (p- mVV ' £ + C\tf\ c °.t> |^| _/9 [(i? _1 |?7|) a/3 (l + |£| P R P ) + exp(-c|£|i?) + |£| a/ W- 
Let us choose R = 2(c|£|) _1 log to obtain 

wt, v (y) -wt <r ,(y) < •l + C , |V’|co.' 3 lfl' 3( “ _1) l ? ?r /3 (iog| 5 p|)- (4.15) 

This implies an estimate for JL — l i as well by evaluating for y G dP, t + \R£. 

~P ~ M < \{V- ~ M) + C\if\ C o,e\Z\K a -%\ af> (iog |^|). 

Here we have used (4.11) to estimate /i — and JL — up. r/ on dP, c + \R£ 1 the error is of the same order 
as in (4.15) so we combined terms. Rearranging the last inequality and making a similar argument for 
the lower bound, we conclude that 

IP - Ml < m c °Mtf {a -%\ a0 (log • (4.16) 

But now we can plug (4.16) back into (4.15) and obtain for any 0 < y • £ < R, 

- w t,ri(y)\ < c'IV , lco.HCI /3( “" 1) l'7r^(iog|y), 


(4.17) 
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the same estimate is obtained for y ■ £ > R by using (4.16) in combination with (4.11). Thus we obtain 
(4.17) for all y £ P$. 

□ 


4.3. Step 3: Reduction to a two-dimensional Problem. The third step of our reduction procedure 
is actually more of notation change. Let F be a homogeneous, uniformly elliptic operator. We are 
concerned with the solution of, 


F{D 2 W£^ rj ) =0 in P ,t 

w i,v = m dy ' v) on dP/: 


(4.18) 


for a fixed unit vector p £ S d 1 with rj ■ £ = 0. 


In the previous section we have already observed that varies only in the £, p directions. To emphasize 
the two-dimensionality of let us define W,£ iV : Ri —> R. by, 


W^ v (z) = W( tV (ziri + z 2 £). (4-19) 

Now W^, n will solve an equation in the upper half space with an operator G v ^ which is essentially the 
projection of F onto the £-77 plane. Let M £ A4 2X 2 a symmetric 2x2 matrix, the definition of G^ tV (M) 
is somewhat cumbersome in terms of notation but the idea is quite simple, 


'■= F ( fj) with fi=V, f2=£- (4-20) 

l<i,j<2 

It is quite important to note the dependence of G on the orientation of ?y; G n ^ may not be the same 
operator as G _ v ^. 


Lemma 4.7. Let Wand G£ :V 

of 


be as given in (4.19) and (4.20). Then W^ v (z) is the unique solution 


Gt-^DlW^) = 0 in 
W^ >v = on <9R)j_. 


(4.21) 


The key point of this reduction is that we realize w^ :V as the solutions of different pdes in the same 
domain with the same boundary conditions. 


4.4. Step 4: The directional limits of y at rational directions. We are now ready to precisely 
characterize the limiting behavior of yf,ip,F) near a rational vector £ £ Z d \ {0} of y(-,ip,F) in terms 
of the boundary layer tails of the class of simpler two dimensional problems (4.21). 

Definition 4.8. Let F be a uniformly elliptic operator as given in Section 2. For if) £ C°^(T d ), a 
rational vector £ £ Z d and a unit vector rj J_ £ define, 

LzivA^’F)) := lim W^r,(0, R; (if), F)). 

R—t 00 

Similar arguments to those used in the previous section will show that L £ is continuous in 77. For example 
see the proof of Proposition 4.3. 

Lemma 4.9. For £ £ Z d \ {0}, 0 £ (0,1) and any 77 , 7 / unit vectors orthogonal to £, 

\L((w, Wh-F)) - (ip,F))\ < C(A, d, /3)\\ip\\ C o,e\p - + log ). 

A combination of the results of the previous sections yields the following classification of the asymptotic 
behavior of y(-,ip,F) near £: 

Proposition 4.10. Let £ £ Z d \ {0} be irreducible and let v : [0,1) — > S^ -1 a geodesic path with unit 
speed and v(0) = £. Then there is a(A) £ (0,1) such that for any f3 £ (0,1), 

\p(v(t),if),F) - L € (i/(0); (t/>,F))| < C(A,dM c o,, (rd) \^t^. 
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5. Continuity of p 

One immediate consequence of Proposition 4.10 is a continuity result analogous to Theorem 4.1 of Choi 
and Kim [16] for operators F(M,y) such that F is rotation invariant. Let us repeat that the proof 
we have given of this result is not new, rather we have made each of the steps [16] quantitative and 
elucidated the secondary two-dimensional cell problem underlying the limiting behavior near rational 
directions. This additional work will be essential to the results that follow but is not so important just 
to get a continuous extension of p{-,ip,F) to the rational directions without an explicit modulus. 

Theorem 5.1. Let £ £ Z d \{0} be irreducible. If F is invariant with respect to the rotations/reflections 
that preserve £ or F is linear, then L^{-\ ( ip,F )) = L^{ip,F) independent of the approach direction. As 
a consequence, p{-,ip,F) extends continuously to £ with value L^{ip,F) and, 

I »M,F) - LsW,F )I < C(A,d,/3)M C o,, (T<i) |£- m a0 , 

for some a(A) £ (0,1) and any j3 £ (0,1). In particular, if F is rotation invariant or linear, p(-,ip,F) 
extends from S d ~ 1 \ WL d to a continuous function on S^ -1 . 

Proof. It suffices to show L ^ is constant in the cases claimed, the rest of the Theorem will then follow 
from Proposition 4.10. For any 771,772 _L £ let O be a rotation sending 771 to 772 and holding £ fixed. Now 
uJ£ i?J i (O*-) has the same boundary data as w^ iV2 (-) and by the rotation invariance of F they solve the 
same pde in P^. Thus by uniqueness they are equal. In particular they have the same boundary layer 
tail so ^( 771 ; = L ? (? 72 ; (tp,F)). 

In the second case we refer to Lemma 3.6 of [19] which shows, using Riesz Representation Theorem, that 
when F is linear and homogeneous p{v,ip,F) = {ip) (the average over the torus). We apply this to 
which satisfies the assumptions of the Lemma since it is a solution of F which is homogeneous and, by 
assumption, linear. We derive for every 77 _L £, 

r i/m 

Ldv, i , <P,F)) = lim w^ v {Ri) = |£| / m ? (t; ( ip,F))dt. 

R-* 00 J 0 

The right hand side is independent of 77 which was the desired result. □ 

As a corollary of Theorem 5.1 we will show an explicit modulus of Holder continuity for the homogenized 
boundary condition when F is rotation invariant or linear. The argument is entirely number theoretic 
and relies on Dirichlet’s Theorem, Theorem 2.11. A sharper estimate in the linear case can be found 
in Section 7. The improvement there is in the rate of convergence at a single rational direction. The 
argument using Dirichlet’s Theorem stays the same. 

Corollary 5.2. Let F satisfying the assumptions of Theorem 5.1. There is a{ A) £ (0,1) such that for 
all fj £ ( 0 ,1), ip £ C°A{ T rf ), and v\ and V 2 irrational vectors in S d ~ l we have 

\p{v u ip,F) - n{u 2 ,ip,F )| < C(d, A, j3)\ip\ c o,f>\u\ - v 2 \ f)a/d . 

Proof. Assume \ip\c°<f> < L the general case follows from scaling. Let e := \v\ — v 2 \, and let N = 
£ -( d ~ l )/ d _ Then due to Lemma 2.11 there exists £ £ h d and n £ Z with 1 < n < N such that 

< (d-l) 1 / 2 ^- 1 ^- 1 ). 

Note that n > | 1 00 1 € I > d _1,/2 |£|. Due to this and the choice of N we have 

\v 2 - n'Vi|oo£| = |ui - u 2 \ + \ v\ - 77 _1 |ui|oo£| 

< £ + (d- 1 ) 1 / 2 n-> 1 | 00 AI- 1 /( rf - 1 ) 

^e + Cdl^N-W-V. 


Vl 


\Vl\ 
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Now we apply Theorem 5.1 with v = Vj at the rational direction £ to conclude that 

1/zOi) - n{v 2 )I < N~ + |^| a / 3 e “/ 3 . 

Using that |£| < N we obtain 

|A*(^i) - ai(i/ 2 )| < Ar#fnog N + (Ae) a/3 < e a/3 / d . 

□ 


6. Discontinuity of h 

Given the set up of the previous sections it may seem at least plausible to the reader that when F is 
nonlinear and not rotation invariant, for a given £ the directional limit function L ^ will typically be 
non-constant, resulting in the discontinuity of the homogenized boundary data. On the other hand it is 
not obvious, at least to the authors, how to prove that any specific pair ( ip, F) results in a non-constant 
Lj. Apart from explicitly computing the solutions the only way to differentiate the boundary layer tails 
of the W^ j7) would be to use maximum principle. However, except in some specially arranged cases, one 
cannot choose ? 7 i, 772 -L £ so that G^ yVl > G^ yV2 and so there is no reason for W^ r] to be ordered in the 
whole of K"} for any such pair 771 , 172 - We instead find monotonicity by perturbing {ip,F). We are then 
able to show that the class of {ip,F) for which L^{-\ (ip,F)) is non-constant is open and dense in the 
appropriate topologies. 

Let us give a heuristic description of how this monotonicity arises. The goal is to show that for any 
( ip,F ) and £ £ Z d \ {0} we can find a nearby {ip',F') such that L^{-; (ip',F')) is non-constant. In this 
paper we are only able to show that a small perturbation of F would lead to L £ being non-constant, 
which directly corresponds to perturbation of the homogeneous operators since F = F. In the general 
case of inhomogeneous F it is not clear to us whether it is possible to perturb F to correspond to the 
desired perturbation of F ; we leave this as an open question. Let us now describe the perturbation of 
homogeneous operators F. First note that we only need to perturb {ip, F ) when L^(-; {ip, F)) is constant, 
otherwise we could take {ip',F r ) = {ip,F). When d > 3 we can find two directions 771,772 perpendicular 
both to each other and to £. We then perturb F in a monotone and hence intrinsically nonlinear way, 
heuristically affecting the choice of diffusions in the 771 direction while leaving the 772 direction unchanged. 
More concretely the perturbation will satisfy that G ^ > G^ ym while G ^ = G^ yV2 . Then, up to a 

small perturbation of ip, strong maximum principle will imply that W^ < W^ yT]1 and, since periodicity 

provides compactness in the lateral directions, also L^{ 771 , {ip', F')) < L^{r] 2 , {ip, F)), while ^( 772 ) remains 
unchanged. Now, having assumed that L^{-; {ip, F)) was originally constant, L{(-; {ip', F')) must be non¬ 
constant. 

The only natural notion of genericity in this setting, to our knowledge, is topological. We make precise 
the topological setting. Our boundary data will be taken from the space, 

C(T d ) = {ip : T d —> M continuous} with the supremum norm. (6-1) 

Let us next define the space of uniformly elliptic operators, 

UEd = {F : A idxd -» K| F £ Ua>i<Sla uniformly elliptic and positively 1-homogeneous}. (6.2) 

Here we recall that Mdxd. is the space of symmetric d x d matrices with real entries. For F £ UE^ we 
define the ellipticity ratio A(F) to be the minimal A > 1 such that F £ <Si,a- It is easy to check from 
this that F £ UE^ are Lipschitz continuous with Lipschitz constant dA{F) with respect to the operator 
norm metric on A id X d- Conversely consider an F which is Lipschitz continuous with respect to the 

d(d+ 1 ) 

operator norm metric on Mdxd — R 5 ■ For this F the gradient DF, from standard inner product 

Tr {AB) on d x d matrices A, B, is defined Lebesgue almost everywhere. The Lipschitz constant of F 
is \\DF\\ L <x>(_M dxd - ) where we implicitly take the underlying matrix norm to be the dual of the operator 
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norm. Based on this definition it is straightforward to check that A (F) < HDFll^. We take as the 
metric on UE^, 

d VEd (Fi,F 2 ) := sup \F,(M) - F 2 (M)\ + WDF, - DF 2 \\ L - {Mdxd y 

||M||=1 

Noting that Cauchy sequences have H-DFnUoo bounded and hence A (F n ) bounded we see that (UE^, duE d ) 
is complete. We draw our operator and boundary data ( ip,F ) from the space, 

X = C(T d ) x UE d with distance d x ((ipi, Fi'), (ip 2 ,F 2 )) = sup \ip! - ip 2 I + d VEd (F 1 , F 2 ), 

which is a complete metric space. 

Theorem 6.1. Let d > 3 and £ £ 1* d \ {0}. Then the set 

= {(ip,F) £ X | y(-,ip,F) is discontinuous at £} is open and dense in X. 

In particular there is a residual set E C X, a countable intersection of open dense sets E = flg e zci\/o}-®£> 
such that for all (ip,F) £ E, 

p(-,ip,F) is discontinuous at every rational direction. 


The proof of the theorem consists of the following two steps. First we prove that E £ is open. The proof 
of Lemma 6.2 is more or less standard, and is due to comparison principle and the stability of viscosity 
solutions with respect to uniform convergence. 

Lemma 6.2. For each t; £ Z d , (if, F) £ X, : {77 £ S d , i) ■ ( = 0} x I -> M is continuous with respect 
to dx at (ip, F), 

sup ( 77 ; (ip',F’)) - Lt.(ry, (ip,F))\ -x 0 as d x ((ip', F'), (ip, F)) -X 0. 

??eS d ,»j £=0 

In particular by Proposition f.lOE ^ is open. 

Next we will show that Ec is dense, whose proof strongly depends on the conditions d > 3 and that F 
is homogeneous. 

Proposition 6.3. Let d > 3 and £ £ Z d . Then for given (ip,F) £ X and e > 0, there exists (ip e ,F e ) 
such that 

d x ((ip e , F e ), (ip, F)) < e and L^(- ; (ip s , F e )) is non-constant. 

In particular y(-,ip e ,F e ) is discontinuous at £ by Proposition f. 10. 


Now we proceed with the proofs. 

Proof of Lemma 6.2. Let (ip n ,F n ) be a sequence in (X,d x ) converging to (ip,F). Let us recall the 
definition of L^(y, (ip n ,F n )) given in Definition 4.8: 


F n (D 2 w n ) = 0 in P£ 

w n (y)=m i (y-y,(ip n ,F n )) on dPf. 


and £{(??; (ip n , F n )) = lim w n (R£). 


(6.3) 


Since F^s are homogeneous, F n = F n but we continue to write F n to emphasize the correct definition 
of L £. We begin by first investigating the continuity properties of me. The claim is 

sup | m^(t; (ip n ,F n )) - m 5 (t; (ip,F))\ -X 0 as n -X 00 . 

t 

Observe that by maximum principle, 

|m 4 (-; ( ip n ,F n )) - m s (-; (ip,F n ))\ < \\ip n - ip\\ oa . 

Thus it remains to show that sup t | m^(t] (ip, F n )) — m^(t; (ip, F))| —x 0. The pointwise convergence with 
fixed t is due to stability of viscosity solutions with respect to uniform convergence of F n , but a little 
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extra work is required to show that the convergence is uniform in t. Note that by Lemma 4.2, we have 
some modulus Q depending on the continuity modulus ui of ip and A (F n ) so that, 

I mt(t; ( ip,F n )) - m £ (t'; {ip,F n ))\ < w(\t - t'\) and {ip,F n ))\ < Halloo- 

Since F n —> F is a convergent sequence in d\jp, d , H-Fnlloo and A (F n ) are bounded. Since m^(-; (ip, F n )) are 
uniformly bounded and equicontinuous j|y-periodic functions on R, every subsequence has a uniformly 
convergent subsequence. It follows that, since m^(-; (ip, F n )) converge pointwise to roc(-;(^,F)), they 
will also converge uniformly. 

Now let us define w n to solve 

f F n (D 2 w n ) =0 in P 5 ; 

\ w n (y) = m^(y ■ r?; (ip,F)) on dP 

Since we have already proven that m^(y-rj] (ip n , F n )) -A m^(y-r]; (ip, F)) uniformly on dP, f, by comparison 
principle Lemma 2.4, 

| lim w n (R£) - 1 ^( 77 ; (ip n ,F n ))\ < sup | w n (y) - w n (y)\ 

fl->oo ye p ( 

< sup \ni^(t] (ip n , F n )) — m^(t; (ip, P))| —)• 0 as n —> oo. 
teR 

By a similar argument as above, since F n —► F uniformly on compact sets of M.dxd when duE d (F n , F ) —► 
0, we have that w n —> w locally uniformly in P ,t and 

| lim w n (R£) - L{(r);(ip,F))\-> 0. 

oo 

Combined with the previous estimate this yields that 

l-C'ffa; (^n,F n )) - L^(rj-, (ip,F))\ -> 0 as n -A oo. 

This shows pointwise convergence of L^(-; ( ip„,F n )). Uniform convergence over all unit vectors 77 _L £ 
will again follow from uniform boundedness and equicontinuity of L^ (see Lemma 4.9). 

□ 


Finally we give the proof of Proposition 6.3. 

Proof of Proposition 6.3. Let £ £ 7L d and ( ip,F ) £ X. If Lc (•; ( ip,F )) is non-constant then we are done, 
so we suppose it is constant and construct (ip e , F e ). 

Let us first show that we can assume without loss that m^(-; (ip,F)) is non-constant. We will choose ip' 
with || ip' — ip\\c(T d ) < ^ such that m^(-; ( ip’,F )) is non-constant. If m^(-; ( ip,F )) is already non-constant 
then we don’t need to do anything and can take ip' = ip. Otherwise let us take 

ip'(y) ■— ip(y) + ecos (2ny ■ £) which satishes \\ip' — ip\\c(T d ) < £• (6-4) 

Observe that for each fixed hyperplane dP,\ + tf we have ip'(y) = ip(y) + £ cos(27r|£|t). Therefore 

m^t; (ip', F)) = mt(t; (ip, F)) + e cos (2 tt |£|£). 

This is evidently non-constant when m^(-;(ip,F)) is constant. If L^(--,(ip',F)) is non-constant we are 
done, thus we may suppose without loss that it is constant. 

Let 771,772 be unit vectors such that rjj _L £ and 771 _L 772 , which is possible since d > 3. We will aim to 
perturb F to construct (F e ,ip e ) so that 

A?(t?i;(< Pe, F e)) < L((rfr-,(ip e ,F s )). (6.5) 


To this end let us define, 


F e (M) := ma x{F(M),F(M + e(ry[Mr]i)r]i <gi 771 )}. 
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It is not difficult to check the definition of ellipticity to see that F e £ UE^. Also, 

sup | F(M + e(r]^ Mr]i)r]i 0 771) — F(M) \ < A (F)e, 

||M||=1 

and furthermore, since D[F(M + s(r]f Mrft)rft 0 771)] = DF + e(rlp DFrft)rh 0 771, 

| DF e — DF | <d eA(F) where it is defined. 
Combining these two estimates it follows that F e is close to F in d\jE d metric, 

dvv d (F,F s )< d A (F)e. 


Let us now show ( 6 . 5 ). From the definition of the 2 -d operators Gj ie for N £ A^2x2, 


G l, V 2 (N) = max |F N t:j a t 10 a-, J ,F Ng(a t g aj + e(m ■ ■ Q:j)t7i 0 771) 

with a\ = 772 and ct2 = £. 

note that 771 • ctj = 0 for j = 1,2 since £, 772 A ??i and so 

Gl V 2 (N) = G^(N) for all e > 0. 

On the other hand, calling e\ = ( 1 , 0 ), 


( 6 . 6 ) 


G l, Vl {N) = max{G i>Vl (N),G^ Vl (N + eNuei 0ei)}. 

Now for any N with Nn < 0 by uniform ellipticity, 

G^ tV 1 (N + eNnei 0 ei) > G^ tVl (N) — eIVu, 

and thus 

G\ m (N) > Gt, m (N ) if N11 < 0 . ( 6 . 7 ) 

We will make use of these observations below. 

We now aim to perturb ip to ip e so that m^(t, (ip e , F e )) = m^(t; (ip, F)). We can always write, 

m^t, (ip,F)) = mt(t, ( ip,F e )) + (m^(t, (ip,F)) - m^(t, (ip,F e ))). 

From the proof of Lemma 6.2 (m^(t, (ip, F)) — m^(t, (ip, F e ))) converges to zero uniformly as e —► 0 . 
Define a modified boundary data ip e , satisfying ip e — ip -A 0 uniformly as e —>• 0 , 

i>e(y) = ip(y) + (m^(y ■ i,(ip,F)) - m^(y ■ £, (ip,F e ))) and we claim m £ (-, (ip e ,F e )) = m ? (-, (ip,F)). 

The claim is immediate from the fact that the added term is constant on hyperplanes parallel to dPc, it 
is the same argument as for the previous perturbation ( 6 . 4 ). Now let us refer to 

Wi r ,(-)=W^(-,(iP e ,F e )) and W^(-) = W £ ,„(-; (iP, F)), 

these two solutions have the same boundary data on dR 2 . but the interior operators are G| and G^ 
respectively. 

From above discussions we know that, 


Wl m < W Lm and Wl V2 = 

where the first inequality is due to the fact G| > Gj i7)1 . On the other hand, since G^ i?)1 (D 2 W| t)i ) < 0 , 
the dichotomy holds due to the strong maximum principle: 

(*) W lvi < or (**) W lm = W £’ r ii • 

In case (i), since and Wp are jiy-periodic in the Z\ direction, 

Wp m (z 1, 1 ) < Wc m (zi, 1) — S for some <5 > 0 . 


( 6 . 8 ) 
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By maximum principle it follows 

Wl m (z) < W^ m (z) - 8 for z 2 > 1 

and therefore 

L?(m; (4>e,F E )) = lim W[ (z) < lim W^ m (z) - 8 < L^m; (ip, F)). (6.9) 

Z2 ~>+00 s ’ Z2 ~>+00 

We have just shown that (6.8) (i) implies (6.9) and therefore we actually have the dichotomy, 

(i) L dm; (tp s ,F e )) < L ? (r?i; (ip,F)) or (ii) W^ Vl (•; (ip e , F e )) = W( tVl (■; {ip, F)). (6.10) 

In case (i) we have that 

mm; (i>e,F e )) < L dm; (i>,F)) = mm; 0 >,f)) = mm; (i>e,F e )) 
which achieves the result since (ip £ , F e ) -A ( , F) in dx as e —> 0. 

We just need to justify that (*) holds. Suppose this is not the case and (ii) holds in (6.10). By Nirenberg’s 
Theorem (Theorem 2.5) II'^. r/l , are C 2 ’ a ° for a small «o(A) > 0 and are classical solutions of their 
respective equations. But then D 2 W^ m = D 2 W^ }m and 

Gi,vi (-)) = 0 = G l m ( D2w l v i 00) = G l m (-° 2 W^ m (z)) for every z G R 2 + 

and therefore, by (6.7), D 2 1 W^ !Vl (z) > 0 for every z € . 

On the other hand, by the 1/|£| periodicity of W^ tV in the z-\ variable, 

ra+i/m 

/ (zi, z 2 )dzi = 0 for all a € K, z 2 > 0. 

J a 

and therefore = 0 in the whole of R+. In particular for any t > 0, W^ Vl (-,t) is constant, 

so W^ V1 is also constant in {z 2 > t} by uniqueness of the bounded solution of G^ :Vl (-) = 0. Since 
t > 0 was arbitrary W^ Vl is constant in R+, but this contradicts the boundary data m^(-; (if>,F)) being 
non-constant. This completes the proof. 


7. Improved Estimates in the Linear Case 


In this section we show the best possible continuity estimates for y(v, ijj, F) by our current methods in 
the linear case. The main tool is the higher regularity estimates available for linear operators in or 
in half-spaces with smooth boundary data. For our purpose W 3,d estimates would be sufficient, but we 
do not pursue this minimal assumption since it would be too much to expect in the general nonlinear 
case anyway. 


Consider u e solving, for v £ S d 1 and R > 1 , 

j F(D 2 u e , |) := —Ty(A(^)D 2 u e ) = 0 in 0 < x ■ v < R 
1 u e = g(x) on x ■ v S {0, R}. 


(7.1) 


We assume that A is Z d -periodic and smooth and satisfies (Id)d X d < A < A(/d) dxd . Due to the linearity, 
the interior corrector can be written as v(y,M) = Y,ijVij(y)Mij, where Vij(y) solves 

-Tr [A(y)(ei ® €j + D 2 Vij)\ = - Aij in M d , 

with the estimate 

KIU ||-D%Hoc <C(A,d). (7.2) 


The following result is not optimal in terms of the required regularity of g , but it is sufficient to improve 
the estimate of Section 4.2, to match the order of the estimate in Section 4.1. 
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Theorem 7.1. Let g and u £ as given above, and let u be as given in (2.4) solving the homogenized 
equation. Then u e converges to u with the convergence rate 

sup \u e (x) — u(x)\ < C(A, d, (3)\\g(R-)\\c3,f>R- 1 £ for any R> 0. 

0 <x-v<R 

In order to use the above theorem to obtain the desired interior homogenization rate, higher regularity 
of the boundary condition is needed. In our setting in Section 4.2, that boundary condition is m^(x ■ fj) 
for some rj _L £. In the following lemma we show that mj(-; (ip,F)) has sufficient regularity to apply 
Theorem 7.1 when if is regular. Indeed the proof of this lemma is the most interesting and delicate part 
of this section. 

Lemma 7.2. Let F be as in (7.1). Then for all k £ N, /3 £ (0,1], 

IK(s < C{k, P, A, d, ||A|| cfc+1 )(log(2 + |e|)) rfc+/31 M w 

Note that by Lemma 4.2 when k = 0, /3 £ (0,1) an improved estimate holds without any logarithmic 
growth in |£|. It is important here that the estimate is not getting too much worse in terms of |£| as k 
increases. 

Now we can return to the proof of Lemma 4.5 and use Theorem 7.1 to achieve the following result: 
Lemma 7.3. Let £ £ Z d \ {0} and g _L £. Define w^, v ,w^ tri as in (4.10). Then, when F is linear, 
sup|wf,„ - W M | < C(/3,A,d, IIAIIcOIKCrMlIca^R^IMtl + ( lo Slfc|| 3 ]- 

p e 

Note that by Lemma 7.2, 

IKKICr 1 ') -w«(0)IIc3./3(r) < f. c0 ,:, < A 11V*lie 7 ) 

where we have used the y|y-periodicity to estimate the lower order terms by \D :i rn^\ C 3..n and we also 
used (log(2 + |£|)) 4 /l£| 3+/3 1$ 1 f° r ]£] > 1- Combining the estimates of Theorem 7.1, Lemma 7.2 and 
Lemma 7.3, in the same way as we did before in the general case in Theorem 5.1 and Corollary 5.2, we 
now obtain the continuity estimate, 

Theorem 7.4. If F is linear andip £ C 7 (T d ), then for every irreducible £ £ Z d \{0} and v £ S ,d_1 \KZ d , 
HFip,F) - L ( (ip,F) | < C(A,d, ||A|| C 5)||V'||c 7 I^ - l£M[l + (log |g_jg H ) 3 ]- 

Furthermore, for every v , v' £ S d ~ x \ RZ d , 

- n(i/,i/;,F) | < C(A,d, ||A|| cs )||V’||c’> - u'| 1/d [l + (log j U ^) 3 ]. 


We omit the proof of Theorem 7.4 as it is a straightforward consequence of the improved estimates of 
Lemma 7.2 and Lemma 7.3 and its usage in the proof of Theorem 5.1 and Corollary 5.2. Before we 
proceed with the proofs of Theorem 7.1 and Lemmas 7.2 - 7.3, we make an interlude to explain some 
background results that we will make use of. 


7.1. Background Results. For the proofs in the next subsection we will need to use the regularity 
results of Avellaneda-Lin [9] for solutions of non-divergence form linear homogenization problems. We 
state the result here in the form that we will use it. Suppose that u s ,v e solve respectively, 

^,,,^2 ,, , . D ^ f -Tr(A(f)D 2 v £ ) = f(x) in B 1 n P v 

-Tr(A(f).DV) =f{x) m B 1 and l 

[ v £ = g{x) on dP„ fl B\. 

The following results hold uniformly in v G S d ~ l . We first state the classical results in unit scale. 
Theorem 7.5. For every fd G (0,1) there exists a constant C{d,A,/d, ||HA|| iD o( T d)) such that, 


(7.3) 


(1) \\D 2 U 1 \\ C o,f l ( Bl/2) < CfosCB! u 1 + ||/||c°.0(Bi)) 

(2) \\D 2 v 1 \\ C o,i3(B 1/2 np„) < C(osc Bl np„ v 1 + ||/||c 0 ^(Binp„) + Il5llc 2 ^(ap^nsi)) 
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Scaling arguments yields that the best possible uniform in e regularity estimate is C 1 ’ 1 . [7,9] shows that 
this estimate indeed holds, although for our purposes the more useful estimate will be the W 2,p estimate. 
Below is a combinations of the main Theorems in [7,9]: 

Theorem 7.6 (Avellaneda-Lin). For every 1 < p < oo, f3 > 0 there exist constants C±(d, A, /3, WDAW^nd^) 
and C 2 (d, A,p, ||Z?A|| iD o( T d)) such that for alls > 0, 

(!) II D 2 u e \\ L ~ (Bl/2) <Ci(osc Bl u e + ||/|| c o,3 (Sl) ) 

(2) ||T> 2 u £ || iP(Bl/2 ) <C 2 (osc Bl u e + ||/|| iP(Bl) ) 

7.2. Proofs of the results of Section 7. Finally we prove Theorem 7.1, and Lemmas 7.2 and 7.3. 

Proof of Theorem 7.1. It suffices to consider the case R = 1, the case for general R > 0 follows from 
rescaling. Let u be as given in the theorem. Let U := {0 < x ■ v < 1}. By the classical C 2,1 estimate up 
to the boundary (for the Laplacian) at unit scale, 

\\D 3 u\\l°°(u) < C(A, d,(3)(oscu+ \\D 3 g\\ c o,p( dU )) < C\\g\\ C 3,p( du y 
Let u e (x) := p e *u(x) with a standard mollifier p e , well defined on U e = {e < x ■ v < 1 — e}, and define 

f> e (x) := u e {x) + 5Z e2 %(f )D 2 ijU e (x). 


Then we have 

W{x) -u{x )| < £\\Du\\ Lao(u) + C(d)£ 2 \\D 2 u e \\ Lao{u) supll^Hioo^d) 

ij 

< C(A, d)\\g\\c20(du)£- 

On x ■ v = e we have, by the regularity up to the boundary of Theorem 7.6, 

\4> e (x) — u e (x)\ < \<f> e {x) — u(x)\ + \u(x) — g(x')\ + \g(x') - u E (x)\ < C{A,d)e\\g\\ c ^o. 

A similar estimate holds on x ■ v = 1 — e. Thus we can estimate, 

sup \u-u e \ < sup \u-u e \ + C(A, d)e||g||c 2 ,i < sup \(f e — u e \ + C(A, d)e\\g\\c^^. (7.4) 

u u e u e 

It remains to estimate sup B \(j) e — u e \. In U e , using the uniform ellipticity and then the definition of the 
corrector, 

-Tr(A(f )D 2 c)>*(x)) > -Tr[A(f )(D 2 tF (x) + £ D 2 Vij ^)D%vT (*))] • • • 

ij 

- Yl A [ 2£ \ DVi i(e)\\ DD ij^(x)\+£ 2 \ v i j (^)\\ D2D i j U% x )\\ 

> —Tr(AD 2 u(x)) — C\D 2 u(x) — D 2 u e (x)\ ■ ■ ■ 

-Ce(sup ||Du ij || 00 )|D 3 u e (a;)| - Ce 2 (sup ||u ij || 00 )]|U 4 ?I e (a;)|. 

Recall that, from (7.2), \\D k v i j\\ 00 < C(A,d) for k = 0,1. Since —Tr (AD 2 u{x)) = 0 we get the following 
supersolution/subsolution conditions, 

f(x) > -Tr (A^)D 2 <j>*(x)) > -f(x) 

where we have called /( x) to be the error from the preceding calculation, 

f(x) := C( A, d)[\D 2 u(x) — D 2 u e (x)\ + e\D 3 u s (x)\ + e 2 |H 4 ?I e (a;)|]. 

The first and third terms can be estimated in terms of u by standard mollification estimates: 

\D A u e (x)\ < e _1 ||iD 3 u||Lcx. and \D 2 u(x) — D 2 u e (x)\ < e\\D 3 u\\l°° for all e<x-v<R — e. 


CONTINUITY AND DISCONTINUITY OF THE BOUNDARY LAYER TAIL 


31 


Thus we have, 


ll/IU- <C(A : d)\\D 3 u\\ L ~s. 


Finally, using C(A, d) ||y ||oo ■ u)(l — x ■ v) as a barrier 

Sup \(f) e -U e | < Sup \(p e - U e | + C(A, d)||/||L=° ( !7 e) , 
u e au e 

which, when combined with the estimate near the boundary gives, 

sup | (f> £ u £ | < C(A,d,p)\\g\\ C 3,f!( du) £. 

U £ 

Combining with (7.4) completes the proof. □ 

The proof of Lemma 7.3 is almost exactly the same as the proof of Lemma 4.5 except we now use the 
improved interior homogenization rate given in Theorem 7.1. 

Proof of Lemma 7.3. Recall that m^(y ■ fj ) is y|y-periodic on dP £ in the direction fj and constant in the 
directions orthogonal to fj. As usual we can restrict to the case |£||??| < 1/2. Let W£ tV , w, y and p as 
given in the proof of Lemma 7.3. Recall from Lemma 3.1 and Lemma 3.2 that, 

\w£,r,(y) ~ Ml + \u>(,v(y) ~~P\ ^ C'[(oscTO { )exp(-c|^|i?) + ||Dw£ >7) || ic » ( p 5) H] for y e dP s + R£. (7.5) 
By Theorem 7.6, for any r > 0, 

||-D^ i?7 || L «, (B|4| _ in p s) < C(|C|oscto 5 + ICrl^wdloo) < q|m £ (|d- 1 -)llc3,,|^| ! 

Fix an R > |£| _1 to be chosen and consider, 

«Wy) = w i>ri (y) + (p-n)R~ 1 y ■ i+ sup [\w itv (-) - y\ + |w £jr) (-) - m|] . 

op s +r£ 

Note that with this modification w^ still solves the same equation as in P^ with the same boundary 
condition on dP^ but now also, 

w^riiy) > wz,r>(y) on dP t + R i- 

Now Theorem 7.1 implies that, 

wt,r,(y) - ^ CWmniR-jWcs.pR^ 1 ^] < C'||m £ (|e|^)||c3./.(i2|^|) 8+ ^- 1 |»7|. 

Rewriting this in terms of iuj 

w t,v(v) ~ w t,v(v) ^ (p-M)-R _1 y-| + C'll^|(l^r 1 ')llc3./3[-R 2+/3 |d 3+ ^l ? ?l+C'(oscm £ )exp(-c|di?) + |^||??|]. 
Let us choose R = 2(c|£|) _1 log to obtain 

U>t,v(y) - w Z,ri(y) < •C + C'||mf(|er 1 -)llc3.ddl ? ll(log]^) 2+/3 - (7.6) 

Due to (7.5), evaluating (7.6) for y £ dP^ + yields 

P ~ P < \{p ~ p) + C'||m ? (|Cr 1 -)llc3^ Idlddog^) 2 ^. 

Rearranging the last inequality and making a similar argument for the lower bound, 

\p~p\ < C'll^(ld _1 -)llc3.d^ll ? 7l( lo glfnF|) 2+/3 - ( 7 - 7 ) 

But now we can plug (7.7) back into (7.6) to arrive at the desired result, 

I wt,v(y)~ w nAy )I < C'l|m ? (icr 1 -)llc3,ddld( 1 0g^) 2+/j - ( 7 -8) 

□ 
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For the proof of Lemma 7.2 we introduce a special norm for functions / : P ed —> R which is suited to 
solving the equation — ]TA A i j(y)D 2 :j u = f(y) in the half-space. See Appendix A for the full develop¬ 
ment. For each R > 0 we define Pp, the class of axis aligned cubes contained in P ed with side length R 
and distance to dP ed at least R/ 2, 

= {Q cube : Q = [-R/2, R/2) d + x with Xd > P}. 


Then we define, 

M p (f, R) = sup (^- f \f(y)\ p dy) and I p (f) = M p (f, 0) + £ N 2 M p (f,N), (7.9) 
q^rWQIJq / 

These norms measure an L p -averaged decay of / in the yd direction. One can easily extend this norm 
to functions / in other half spaces by simply rotating appropriately the cubes used in the definition. 
We will use the norms M p (f 1 R) and I p (f) for functions / : P„ — > R for a general v £ 5' d—1 with¬ 
out further comment. Roughly speaking, when I p (f) is finite then there exists a bounded solution of 
— Aij(y)D?jV, = f{y) in P„ with zero (or a bounded if) Dirichlet boundary data and one can bound 
supp |it| < Again we refer to Appendix A for the full explanation. 

Proof of Lemma 7. 2. Recall that ( t) is defined as the boundary layer tail of td solving the cell 

problem, 


-E ij A iAy + T ) D ij v Lr = o 
h, t = ip( y + T ) 


on 


dP f 


(7.10) 


Let us denote the differential operator d := £ ■ D T . Take any p > d/ 2, for example p = d will work fine. 
We will prove by induction that, for all k £ Nil {0}: if ||^||c fc + 3 < 1 then the following hold for all R > 0 
uniformly in r £ R d , 


sup \d k v^ T (y)\ < k , A (log(2 + \f\)) k 
y&Ps 

osc d k v/:, T (y) < kA (log(2 + \f\)) k exp(-c 0 P/|£|) 
yd>R 

M p (D 2 d k v^ T ,R) < k , A (log(2+ |£|)) fc (2 + P)- 2 exp(-c 0 P/|f|) 
We note that if (in) holds for some k then by Lemma A.3, 

I p (D 2 d\ tT )< k> A (l°g(2 + |£|)) fc+1 


W 

(**) 

(in) 

(iv) 


Once we have proven that (i) holds for all k we will be done since = limp-^ d k v^ t g(R£). The 

outline of the argument is as follows. For each k we prove ( i)-(ii ) using that ( i)-(iii ) hold for all m < k. 
Then we show (Hi) using Theorem 7.6 to estimate M p (f , R) and thereby I p (f)- 

For k = 0, (i) is maximum principle and (ii) is a consequence of Lemma 3.1. To show (Hi) note that 
Theorem 7.6 implies 


\D\,r\ 


L°°({yt>R}) 


<A R 


-2 


OSC V£ t 

{yi>R/ 2} 


<a(R + 2) 2 (oscf>) exp(— ic 0 i?/|£|) for R > 1, 


where Co is the exponential rate of convergence to the boundary layer tail from Lemma A.4. For R < 1 
the estimate follows from the up to the boundary C 2 ” 0 estimates at unit scale, see Theorem 7.5. 

Suppose that (i)-(iii) hold for all m < k — 1. We aim to prove (i)-(iii) for d k v^ tT under the assumption 
C k+3 < 1. Our induction is based on the fact that d k v^ tT solves 

- 5Z a h(v + r) D ld k v LT (y) = kAv) 


in Pjt, 


(7.11) 


CONTINUITY AND DISCONTINUITY OF THE BOUNDARY LAYER TAIL 


33 


where 

kAy)~Yl [(£• D v) eA ij(y + T)\ D i j d rn v itT {y) 

ij £-\-m—k,iy^0 


with boundary data 

d k ViAy) = [(£ • D y )A\(y + r) on 8P,. 

(To be completely precise we should take difference quotients instead of 9 fc uj jT ; the reader can easily see 
how to make our formal argument rigorous.) Note that the boundary data d k v^^ T , and the operator 
F are all periodic with respect to a lattice on dP,\ with unit cell diameter of order |f|. This will allow us 
to use the results of the Appendix which extend Section 3 Lemma 3.1 to include equations with a right 
hand side. We use the inductive hypothesis (in) to see that f$Av) satisfies, 


M p (At, R) <k,A (log(2 + |f|)) fc 1 exp(—coi?/|f |). 


In particular f^ T fits under the assumptions of the results of the Appendix. We can apply directly 
Lemma A.3 in combination with Lemma A.2 to get, 

sup d k v^r(y) < ky A II3VIU + Ip(h,r) < (log(2 + |f|)) fc . (7.12) 

yeP 4 

We have used that, by assumption, ||^||c fc + 3 < 1- Then we use Lemma A.4 to get, 
osc d k v itT <k,A [Halloo + (log(2 + |f|)) fc ] exp(-co-R/lfI) 

y-v>R 

< (log(2 + |f|)) fe exp(—co-R/lf |). (7.13) 

This establishes (ii). 

It remains to prove ( Hi ). Theorem 7.6 yields that, for R > 1, 

M p (D 2 d k v i r,R)< R~ 2 osc vt r + M p (f J T , §) 

{vi>R/ 2 } 

< k ,A (2 + R)~ 2 ( log(2 + |f |)) fe exp(—c 0 i?/|f |). 


For R < 1 the estimate follows from the up to the boundary C 2 ^ estimates at unit scale, see Theorem 7.5, 
combined with (7.12) and ||'0||c fc + 3 < 1- Together we obtain for all R > 0, 

M p (D 2 d k v^ T ,R) < k ,A (2 + I?) _2 (log(2 + |f|)) fe exp(—co-R/|f|). (7.14) 

□ 


Appendix A. 

Here we prove Lemma 3.1 but now including a right hand side in the equation. We restrict to linear 
equations, although a version of these results for nonlinear equations is true as well with a stronger 
assumption on the right hand side. Let us consider v solving the following problem: 

( -T,ij A ij(y) D ijV = f(y) in p ed 
1 v = <f>(y) on dP ed . 

We assume that A, f and <j> are periodic with respect to linearly independent translations li, ...id -1 £ 
dP ed . Recall that we denote Z by the periodicity lattice generated by {£j }, Q its unit cell, and L the 
diameter of Q. Let us assume L > 1. We will furthermore assume that A is periodic with respect to 
some rotation of the Z d lattice and A £ C°’ /3 (K d ) for some /3 > 0. This periodicity and regularity will 
only be used apply the results of [9] to obtain estimates for the Green’s function G(x, y) associated with 
the inhomogeneous operator and the domain P ed . 

When / is continuous and has sufficient decay there is a unique bounded solution of (A.l). In order to 
investigate v we define an appropriate norm for the decay of / at infinity in the yd direction, this norm 
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measures the L p -averaged decay of / far from dP ed . For each R > 0 we define Tr, the class of axis 
aligned cubes contained in P ed with side length R and distance to dP ed at least R/ 2, 

J-r = {Q cube : Q = [—72/2, R/2) d + x with Xd > R}. (A.2) 

Then we define, 

M p (f,R) = sup (-3- [ \f(y)\ p dy) ** and I p (f) = M p (f, 0) + £ N 2 M p (f,N), (A.3) 

Q^rWQIJq J ^ 

where we define M p (f, 0) = ||/||l°°(,p )• It will also be useful to define I p (f,R) where the sum in the 
definition of I p (f) is restricted to 2 N 3 N > R, 

I P (f,R) = E N 2 M p (f,N). (A.4) 

2 ®3N>R 


Now we remind about some facts about the Green’s function for operators with periodically oscillating 
coefficients. Recall that for each x £ P ed the Green’s function G(x, y) for the domain P ed solves the 
adjoint equation, 

{ -Y /i j D yi D y j ( A ij(y) G ( x ,y)) = S (.y-x) for y£Pe d 
\ G{x,y) =0 for y £ dP £d . 

From the work of Avellaneda-Lin [7-9] G{x, y) satisfies many of the same estimates as the Green’s 
function associated with the Laplace operator. We will make use of the following, 


Theorem A.l. (Avellaneda-Lin [7 -9]) Let (3 £ (0,1) and suppose that 1 < A{y) < A is periodic on 
with respect to (or a rotation of the integer lattice) translations and ||A|| C o,/j( R <i) < oo. There is 
c = c(d, A, ||A|| C o,/ 3 , /?) < oo so that the Green’s function for the half-space G(x , y) solving (A.5) satisfies, 


G(x, y) < c XdVd , d and G(x,y) < 
\x - y\ a 


c| log \x- ?/|| d = 2 
c\x — y\ 2 ~ d d > 3. 


(A.6) 


Note that the result proven in [8] is actually for divergence form operators but a standard trick, which 
is explained in [7,9] , allows one to translate the results for non-divergence form equations as well. The 
point is that non-divergence form equations can be written as a divergence form equation with uniformly 
elliptic matrix A' satisfying divA' = 0. 

Now we are able to state and prove, with the aid of Theorem A.l, our result on the upper and lower 
bounds of solution of (A.l) in half-spaces. 


Lemma A.2. Let p > | and suppose I p (f) is finite. There exists C{p,d, || A|| c °,^, (3) so that there exists 
a unique bounded solution v of (A.l) with, 

sup \v\ < max \<j>\ + CI p (f) when d > 3 
P e d dP *t 

sup \v\ < max \(j>\ + C[I p (f) + sup N 2 (logN)M p (f,N)\ when d = 2 
Pe d dP *d Ne 2 N 


Proof. It suffices to show the result with f = 0. We just need to show that 



G{x, y)\f{y)\ dy < I p (f) 


(A.7) 


so then 


v(x)= f G(x,y)f(y) dy 
Jp d 

is well defined, solves the desired equation in P ed and has the upper bound claimed by the Lemma. 


Fix an x = (x r , xd) £ Pe d , we partition P £d by dyadic cubes which are adapted to the location of x. This 
is a technical step, the idea is to make the computations later easier by keeping all of the cubes in the 
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partition but one a controlled distance from the singularity at x of G(x,y). If x d > 1 let N x £ 2 N so 
that N x < Xd < 2 N x and call a £ [1,2) such that aN x = Xd , if Xd < 1 then set a = 2. Then we define 
the cubes Qnj = x' + aN(j , 1) + [— aN/2, aN/2) d , N £ 2 N and j £ Z d_1 , with side length comparable 
(by a factor of 2) to their distance to the boundary. Note that by our set up x £ Qn x ,o and the distance 
of x to any other cube Qnj is at least c d N{ 1 + \j\). 

Then we estimate v by, 

[ G(x,y)\f(y)\ dy = f G{x,y)\f(y)\ dy 


<EM 1/P (/ G(x, y y'dy 
Q Q 


1/p' 


\Q\ 


\ i /p 

\f{y)\ p dy I 


S E E N d / p M p (f, N ) 

We claim that for any p > | and every N £ 2 H , j £ Z d_1 : 

/ \ i/p' 


\ i/p' 

G{x,y) p ’dy j 


IQnjI Jq nj 


G(x, y) p dy < C(p, d, ||A|| C o,,,/3)iV^(l + 


(A.8) 


Taking for granted (A.8) we continue the computation 

i / p ' 


E E N d/p M p (f, N )( f G(x, y) p ' dy\ <E E N d ^ + ^ Mp (f, N)N 2 ~ d (l + \j\)~ d 

at e 2 K jez'i- 1 J jve2 N je i d ~ x 


= E E ^ 2 Mp(/,A0(i + bi) 

jVe2 N j£ z^- 1 

<d E N 2 M p (f,N) = I p (f). 

AT£2 n 


-d 


This proves the desired estimate (A. 7). 


To finish the proof we need to justify (A.8). Let y £ Qtvj with either j y 0 or TV y N x . If IV > A^ x then 
|a; — y| > (1 + |j|)JV and so using the first part of (A.6), 

G(x,y)<N 2 - d (l + \j\)~ d . 

If N < N x then Xd~ y d > N x /2 = a;<j/2 and \y' — a/| > A7|j |/2 so using the first part of (A.6) again, 

Nx d 


G(x,y) 


< 


N 2 d \j\ d N\j\>x d <N 2 -d (1 + \ -o-d 
(x d + N\j\y - \ Nx 1 ^ N\j\ < x d ~ ( M) ' 


< 


Finally, if j = 0 and IV = N x then we use the second part of the Green’s function estimate (A.6) and 
p’ < ^2 to g e G when d > 3, 


' Qa 


\ 1 /p' 

G{x, y) p 'dy ] < 


(f \x-y\^ p 'dy] <([ 
\jQ n x ,o ) \Jo 


r VdN \ 1 / p 

r (2-d)p' r d-l^ r j < N 2-d N d/p' 




which, after dividing through by N d / p ~ IQat.oI 1 ^ j is the desired estimate (A.8). When d = 2 we just 
use the d = 2 estimate from (A.6) instead. 

□ 
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For notational simplicity we will only focus on the particular case to be used in the paper: let us assume 

M p (f, R) < K( 2 + R)~ 2 e~ bR with b > 0. (A.9) 

With above assumption we estimate I p (f) in the following Lemma: 

Lemma A.3. Let f such that M p (f,R) satisfies (A.9) for all R > 0, then 

I P (f)<CK log(2+i) and I p (f,R)< CK\og(2+\)e~ bR . (A.10) 

Also, for the d = 2 case, we have, 

N 2 (\ogN)M p (f,N)<CK\og(2+l) and sup A 2 (log N)M p (f, N) < CK\og{2 + \)e~ bR . (A.ll) 

N>R 

Proof. Without loss of generality we can set I\ = 1. Note that by taking R —>■ 0 in (A.9) (with K = 1) 
we see M p (/,0) = \\f\\ L °°(p ed ) < 1. 

w) < 1 + E n 2 m vU,n) < i + e e ~ bN 

NG 2 k JVG2 n 

<log(2+i)+ ^ e~ bN 
N> i 

<106(2 + 1)+^^ 

N>i 

<log(2+i). 

where we have used e x > x in the second inequality. A similar argument proves the estimate of I p (f, R) 
and (A.ll) is straightforward from calculus. □ 


Now we are able to prove the corresponding version of Lemma 3.1 when / is of the form in (A.9). 

Lemma A.4. Suppose that f satisfies (A.9) for some K,b > 0. Then there exists p,(<f),F,f) and 
co(A,d) > 0 such that, 

sup \v(y) - n\ < C(p,d, ||A|| C o,3,/3)[(osc</))e _Coii/i + if log (2 + \)e~ bR/L }. 
ye d >R 

Proof. We just argue when d > 3, but it will be clear given the matching estimates obtained from 
Lemma A.2 and Lemma A.3 that the same result holds in d = 2. 

After rescaling we may assume that K,osc(f> < 1. Let a{d, A) £ (0,1) and Co as given in Lemma 3.1. 
Next define A(j) := 2I p (f,jr) + r 2 M p (f, ( j + 1/2 )r) for j € N and r > 2C^j a L. We claim that 


k -1 


osc V < Co(2L/r) ak + ^ C*- j (2L/r) a ^A{j) + 2I p {f, kr). 
d j =0 

Let us assume for now that (A. 12) is correct and continue with the rest of the proof. Define 

r := (2eC 0 ) 1/ “L and c 0 := (2eC 0 )~ 1/a 


(A.12) 


so that C^{2L/r) ak = e~ k and 


k -1 


osc v<e k +e e J A(j) + 2I p (f,kr). 

P e +krv ~ ^ / 

d 3 =0 

Due to our assumption (A.9) we have the bound 

r 2 M p (f , (j + l/2)r) < r 2 (2 + ( j + l/2)r)- 2 e - b ^ +1 ^ r ' L 

< e ~bjr/L < e ~CQ 1 bj 


(A. 13) 
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and using Lemma A.3 we also have, 

W,jr) < b log(2 + l)e~ b ^ L ~ log(2 + l)e~ c ^. 

Plugging these estimates into (A. 13) yields 

k -1 

osc v < e~ k + log(2 + l)e~ k V e^" 1 ^ + log ( 2 + i) e ~ c o lbk 

Pe^+kris 

d i=o 

< e~ k +log(2+ l)e~ c ° lbk . 

Now for any R > 0, let k = [R/r] to get the desired result 

osc v < osc v < e ~ c ° R / L + log(2 + \)e~ bR l L . 

Pe d +Rv Pe d +krv ° 

It remains to prove (A.12) by induction. For k = 0 (A.12) follows from Lemma A.2. Assuming (A.12) 
for k, we prove it for k + 1. Note that 

Vk(y) = v(y + kre d ) 

satisfies — A^{y + kre d )D 2 ^Vk = f(y + kre d ) in P ed with boundary data <fik{y) = v(y + ke d ). Both 
the operator (by assumption) and the boundary data (by ^-periodicity of operator and uniqueness) are 
periodic with respect to (lj)jZi translations, and 

fc-i 

oscvfe < C k (2L/r) ak + £ C k ~ j (2L/r)^ k ~»A{j) + 2J p (/, kr) 

3=0 

by the inductive hypothesis. Then by the interior Holder estimate in B r / 2 {re d ), 

\vk\c-{B r/i {re d )) < 2 a C 0 r~ a (osc v k + r 2 M p (f , (k + l/2)r)), 
and so, since r > 4L, the oscillation on Q can be estimated by 

osc v k < C 0 {2L/r) a {oscv k + r 2 M p (/, (k + l/2)r)), 

Q+re d 

On the other hand v k is periodic on dP ed + re d with respect to the translations and periodicity 

cell Q so osc Qp Cd+red v k = osc Q +red v k - Using the inductive hypothesis the right hand side above is 
bounded by, 

fc-i 

osc v k < C 0 (2L/r) a (C k (2L/r) ak + V C k ~ j (2L/r)^ k ~AA(j) + 2 /„(/, kr) + r 2 M f ((k + l/2)r)) 

9P ^ +red 

k -1 

= U 0 fc (2L/r)“ (fc+1) +J2 C o +1 ~ J ^L/r) a{k+l -AA{j) + C 0 (2L/r) a A(k) 
j=o 
k 

= C k (2L/r) a(k+1) + J2 C k+1 - J (2L/r) aik+1 - j) A(j). 

j=o 

Finally using Lemma A.2 in place of maximum principle to bound the oscillation in the entire half space, 
osc v = osc v k 

Pe d +(k+l)re d Pe d +re d 

< OSC v k + 2I p (f 1 (k + l)r) 

dPe d +re d 

k 

< C k (2L/r) a ^ + J2 C k+1 ~ j {2L/r)^ k+1 -AA{j) + 2 I p (f, (k + l)r), 

3=0 


which is the inductive hypothesis for k + 1. 


□ 
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